arXiv:1508.04565vl [astro-ph.GA] 19 Aug 2015 


Astronomy & Astrophysics manuscript no. ngc3227 ©ESO 2016 

August 16, 2016 


A Variable-Density Absorption Event in NGC 3227 mapped with 

Suzaku and Swift 

T. Beucherinil, A.G. Markowit2lli^LH F. KraulCI 3 , G. MiniuttP, A.L. LonginottP, M. GuainazzP, I. de La Calle Pere2P, 
M. MalkaiP, M. ElviP, T. MiyajfiniH D. HiriarP, J.M. LopeP, I. AgudcP, T. DauseP, J. GarciP, 

A. KreikenbohnJnH, M. KadleP, and J. WilmP 

* Dr. Remeis-Sternwarte & Erlangen Centre for Astroparticle Physics, Universitat Erlangen-NUmberg, Stemwartstrasse 7, 96049 
Bamberg, Germany 

^ Lehrstuhl fur Astronomie, Universitat Wurzburg, Emil-Fischer-StraBe 31, 97074, WUrzburg, Germany 

^ Center for Astrophysics and Space Sciences, University of California, San Diego, 9500 Gilman Dr., La Jolla, CA 92093-0424, 

USA 

Alexander von Humboldt Eellow 

^ Centro de Astrobiologia (CSIC-INTA), Dep. de Astrofisica, European Space Astronomy Centre, PO Box 78, Villanueva de la 
Canada, 28691 Madrid, Spain 

® Instituto de Astronomfa, Universidad Nacional Autonoma de Mexico (UNAM), 04510 Ciudad de Mexico, D.F. Mexico 
’ European Space Agency, European Space Astronomy Centre, PO Box 78, Villanueva de la Canada, 28691 Madrid, Spain 

* Physics and Astronomy Department, UCLA, Los Angeles, CA 90095-1562, USA 
® Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 

Instituto de Astronomfa, Universidad Nacional Autonoma de Mexico, Km 103, Carret, Tijuana-Ensenada, Ensenada 22860, BC, 
Mexico (P.O. Box 439027, San Ysidro, CA, 92143, USA) 

" Instituto de Astrofisica de Andalucfa (CSIC), Apartado 3004, E-18080 Granada, Spain 

xx-xx-xx / xx-xx-xx 


ABSTRACT 

Context. The morphology of the circumnuclear gas accreting onto supermassive black holes in Seyfert galaxies remains a topic 
of much debate. As the innermost regions of Active Galactic Nuclei (AGN) are spatially unresolved. X-ray spectroscopy, and in 
particular line-of-sight absorption variability, is a key diagnostic to map out the distribution of gas. 

Aims. Observations of variable X-ray absorption in multiple Seyferts and over a wide range of timescales indicate the presence of 
clumps/clouds of gas within the circumnuclear material. Eclipse events by clumps transiting the line of sight allow us to explore the 
properties of the clumps over a wide range of radial distances from the optical/UV Broad Line Region (BLR) to beyond the dust 
sublimation radius. Time-resolved absorption events have been extremely rare so far, but suggest a range of density profiles across 
Seyferts. We resolve a weeks-long absorption event in the Seyfert NGC 3227. 

Methods. We examine six Suzaku and twelve Swift observations from a 2008 campaign spanning 5 weeks. We use a model account¬ 
ing for the complex spectral interplay of three differently-ionized absorbers. We perform time-resolved spectroscopy to discern the 
absorption variability behavior. We also examine the IR-to-X-ray spectral energy distribution (SED) to test for reddening by dust. 
Resuits. The 2008 absorption event is due to moderately-ionized (log^ ~ 1.2-1.4) gas covering 90% of the line of sight. We resolve 
the density profile to be highly irregular, in contrast to a previous symmetric and centrally-peaked event mapped with RXTE in the 
same object. The UV data do not show significant reddening, suggesting that the cloud is dust-free. 

Conciusions. The 2008 campaign has revealed a transit by a filamentary, moderately-ionized cloud of variable density that is likely 
located in the BLR, and possibly part of a disk wind. 

Key words, galaxies: active - galaxies: individual (NGC 3227) - galaxies: Seyferts - X-rays: galaxies 


1. Introduction 

X-ray spectroscopy allows us to probe circumnuclear matter ob¬ 
scuring the direct line of sight to the supermassive black holes 
(SMBHs) that power AGN. This material can be part of a dusty 
torus surrounding the nucleus ( |Antonucm 19931, of the Broad 
Line Region (BLR), or a disk outflow. Variable absorption has 
been found for a number of AGN, e.g., a sample of 12 Seyfert 1.5 
galaxies with inclination angles comparable to the torus open¬ 
ing angle, which included NGC 3227 ( |Beuchert et al.| 2013| l as 
well as other AGN with intermediate optical classifications such 
as NGC 4051 dGuainazzi et al.||1998|l, MCG-6-30-15 ( |McKer- 
|nan & Yaqoob|1998| l, or NGC 3516 ( Turner et al.|2008] ). Risaliti 


et al. ( 2002[ ) study an X-ray selected set of 25 Seyfert 2 galax¬ 
ies and And soft X-ray variability on both long and short time- 
scales. Short-term (~1 d) absorption events have been detected, 
e.g., for Mrk 766 (|Risaliti et al.||2011|l, NGC 5506 (|Markowitz 
et al.|2014|i, NGC 4388 |Elvis et al.|2004l i and NGC 1365 ( |Risal- 
iti et al. 2007 2009a|b|. Long-term events (> 7 d) have been 
found for, e.g., CenA (Rivers et al.||2011ai, Eairall 9 (|Lohfink 


[eraL]|20T2l l, or NGC 3227 in 2000/2001 ( |Lamer et al.pOOgj h 
thanks to the Rossi X-ray Timing Explorer (RXTE) monitoring. 

These absorption events can be explained by transits of dis¬ 
crete clouds or clumps of gas across the line of sight to the 
central X-ray continuum source (e.g., Risaliti et al. 2002[ Lo 


|hflnk et ar||2012[ [Markowitz et al.||20i4[ ). They support a new 
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generation of “clumpy torus” models ( |Elitzur||2007 INenkova] 
et al.|2002 2008a|b I. In these models, clouds are typically con¬ 
centrated towards the equatorial plane hut with a soft-edge an¬ 
gular distribution, usually on near-Keplerian orbits, and possi¬ 


bly embedded in a tenuous inter-cloud medium (Stalevski et al. 
2012|l. Outside the dust sublimation radius, which is typically 


several light-weeks away from the SMBH, the presence of such 
a clumpy component is supported by fits to the infrared spectral 


shape (As ensio Ramos & Ramos Almeida|2009[|Alonso-Herra^ 
jet al.|20lT| l. Closer to the SMBH the population of dusty clouds 

may transition to the dust-free clouds that comprise the BLR 
( jElitzur 2007| l. A clumpy X-ray absorbing medium located at 
distances commensurate with the BLR - and possibly identified 
as BLR clouds themselves - has also been suggested (|Risahti 


et al.|2009a 201 l|l. Arav et al. s 1998 investigation of optical line 


mapping, ~ 20 Id (Suganuma et al. 2006) 1, these clouds are likely 
in the inner dusty torus, or at least the outermost BLR. 

In this paper we present an additional, intriguing, absorp¬ 
tion event in NGC 3227 found in quasi-simultaneous Suzaku and 
Swift data from 2008. The column density profile turns out to be 
highly irregular. In Sect. we present the data reduction while 
Sect, [^revisits the two already published XMM-Newton obser¬ 
vations in order to define a model that we use to simultaneously 
fit the archival Suzaku observations. The X-ray and UV data pro¬ 
vided by Swift are also analyzed. In Sect. we discuss the prop¬ 
erties and location of the absorbing cloud and its possible origin 
in light of these observations. 

2. Observations and Data Reduction 


prohies in NGC 4151 also supports such a distribution of discrete 2.1. Observations 


BLR clouds, with a large number of rather small clouds. Erom an 
observational point of view, however, the mid-IR emission only 
probes matter outside the dust sublimation zone, while X-ray 
absorption probes the full radial range. Many of the short-term 
X-ray absorption events have been interpreted as BLR clouds 
( |Risahti et al.||2007[ ), while the longer term events found, e.g., 
by Markowitz et al.| ( 2014[ durations of more than a few days to 
more than a year) were inferred to be due to clouds residing in 
the outer BLR or the inner dusty torus. 

When adequate data have been available, i.e., sustained sam¬ 
pling on timescales longer than the eclipse duration that re¬ 
solves the eclipses, one has been able to pinpoint their ingress 
and egress, and also to obtain time-resolved information on 
the column density prohle NniO along the transverse direc¬ 
tion. Only a few density profiles have been resolved; The events 
in NGC 3227 and Cen A mentioned above featured symmet¬ 
ric, non-uniform, and centrally-peaked column density profiles. 
An event in Mkn 348 may also fall into this category (|Aky- 


Turner 


tures indicate that some warm absorbers are outflowing along 
the line of sight ( |Blustin et al.|2005|l and m ay be launched from 
the inner accretion disk (Krongold et al. 2007| ). The observed 
warm absorbers could be part of stratified outflows ( |Tombesi] 
jet al.|[20T3] l. Such disk winds can be launched by magneto 
hydrodynamical (MHD) forces ( Blandford & Payne) 1982) Con 


[topoulos & Lovelace) 1 994) [Eukumura et al.|2010|l. 

The Seyfert 1.5 AGN NGC 3227 was subjected to sustained 


monitoring with RXTE from 1999-Jan-02 to 2005-Dec-04 (Ut 


|tley & McHardy||2005| ). The monitoring revealed two eclipse 
events: an ~80-day even t in 2000/2001 (|Lamer et al.|2003 1 and a 
2-7 day event in 2002 ( Markowitz et al.||2014| l. The estimated 
distances from the central engine are tens of light-days (Id). 
Given the estimated location of the BLR from the SMBH, at 
2-19 Id ( |Peterson et al.|[2()04) |Landt et al.||2008] l, and the in¬ 
ner edge of the IR-emitting torus as determined by reverberation 


Eigure[T] provides an overview over all archival observations of 
NGC 3227 that we consider. NGC 3227 was the subject of a 
35 d long Suzaku and Swift monitoring campaign between 2008 
October 2 and 2008 December 03. The campaign consists of 
six Suzaku and twelve nearly simultaneous Swift observations 
(Table [T]|. Outside of this campaign, further archival data with 
good soft X-ray coverage are provided by two XMM-Newton ob¬ 
servations, one from 2000 November 28-29 at the onset of the 
2000/2001 absorption event (obsid; 0101040301, exposure time 
post screening 32ks; [Lamer et al.|2003 1 and a 97 ks observation 
on 2006 December 3-4 (obsid: 0400270101; [Markowitz et al.| 
2009)1 during a state with intrinsically high flux and relatively 


las et al. 2002 1 . Maiolino et al.[ (2010i report comet-shaped 
clouds, with dense “heads” and less-dense “tails”, in NGC 1365. 
)Markowitz et al.)(2014) 1 report a doubly-peaked absorption event 
in NGC 3783. These results indicate a broad variety of profile 
shapes, hinting at a range in cloud origins and/or the physical 
mechanisms that shape clouds. 

To add complexity to the interpretation of X-ray spectra, a 
significant number of Seyferts also show evidence for a com¬ 
plex interplay of differently ionized layers of warm absorbing 
gas along the line of sight (e.g., Blustin et al. 2005) 


low absorption. The observation from 2006 yielded reflection 
grating spectrometer (RGS) data of sufhcient signal-to-noise for 
a detailed analysis of the absorbers. Recent Swift pointings from 
2013 October 20 until 2015 May 02 are also included in our 
analysis. 

Eor the reduction of these data we used HEASOET v. 6.15.1 
and XMM-Newton SAS v. 13.5.0. The spectral analysis is per¬ 
formed with the Interactive Spectral Interpretation System (ISIS; 
Houck & Denicola|2000l. Uncertainties on spectral fit parame¬ 


ters correspond to the 90% confidence level for one interesting 
parameter. Throughout the paper, we use the systemic redshift of 
z = 0.00386 ( de Vaucouleurs et al.|1991| and a G alactic equiva¬ 
lent h ydrogen column of Ah = 1.99 x 10^° cm“^ (Kalberla et al. 


2005 I. The black hole mass is assumed to be 1.75 x 10' Mq, the 


[et al.|2008)|. This gas likely originates in the accretion disk, the 
BLR, or the inner torus ([George et al.[[199^ [Krolik & Kriss) 
2001) [Reynolds & Eabianj 1995) 1. Blue-shifted absorption fea- 


average of various measurements collected at the AGN Black 
Hole Mass Databas^ Luminosities are calculated using a lumi¬ 
nosity distance of 20.3 Mp c ([Mould et al. [2000)1. W e perform the 
X-correction according to Ghisellini et al. p()09) l. Neutral ab¬ 
sorption is modeled with tbnew, an improved version of the ab¬ 
sorption model of Wilms et al.[ ( [2000j l, using cross sections from 
[Verner et al.[ ( [ 1996) 1 and abundances from [Wilms et al.[ ( [2000[ l. 
Ionized absorption is modeled using the zxipcf model, which 
is based on XSTAR tables ( Kallman & Bautista[[2001[ l. In this 
model, following Tarter et al. ( 1969) 1 the ionization parameter, 
for a medium at distance R from a source of luminosity Lion is 
defined as 

f^U,J{nnR^) ( 1 ) 

where hh is the hydrogen number density of the absorber. 

2.2. XMM-Newton 

Both XMM-Newton observations were taken in full-frame mode 
of the EPIC-pn camera ( Striider et al.[200T Turner et al. [2001) 1. 


http://WWW.astro.gsu.edu/AGNmass/ 
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Table 1. Suzaku and Swift observations in 2008 with their screened exposure times. Counts given are in the energy bands given in Sect. The 
check-symbol (y/) denotes observations that are used in the data analysis, while observations labeled with a cross (X) are excluded (see Sect. |2.4| 
for details). 


Suzaku XIS 






Swift XRT 






abbrv. 

obsid 

det 

time 

exp [ks] 

cnts [xIO**] 


abbrv. 

obsid 

mode 

date 

exp [ks] 

cnts 


Suz 1 

703022010 

XISO 

2008-10-28 

58.9 

7.8 

V 

Sw la 

00037586001 

pc 

2008-10-28 

1.0 

543 

V 



XIS 1 



8.6 

V 

Sw lb 

00037586002 

pc 

2008-10-29 

2.0 

2386 

X 



XIS 3 



8.0 

V 










HXD 


47.9 

3.4 

V 








Suz 2 

703022020 

XIS 0 

2008-11-04 

53.7 

2.2 

V 

Sw 2a 

00031280001 

Wt 

2008-11-04 

1.0 

255 

V 



XIS 1 



2.2 

V 

Sw 2b 

00031280002 

Wt 

2008-11-05 

2.2 

497 

V 



XIS 3 



2.2 

V 










HXD 


46.4 

3.1 

V 








Suz 3 

703022030 

XISO 

2008-11-12 

56.6 

3.3 

V 

Sw 3a 

00031280003 

Wt 

2008-11-12 

1.0 

291 

V 



XIS 1 



3.4 

V 

Sw 3b 

00031280004 

pc 

2008-11-13 

2.2 

702 

V 



XIS 3 



3.4 

V 










HXD 


46.6 

3.2 

V 








Suz 4 

703022040 

XISO 

2008-11-20 

64.6 

1.5 

V 

Sw 4a 

00031280005 

pc 

2008-11-21 

2.1 

280 

V 



XIS 1 



1.8 

V 

Sw 4b 

00031280006 

pc 

2008-11-22 

1.9 

256 

V 



XIS 3 



1.5 

V 










HXD 


43.4 

2.9 

V 








Suz 5 

703022050 

XIS 0 

2008-11-27 

79.4 

3.7 

V 

Sw 5a 

00031280007 

pc 

2008-11-25 

2.0 

429 

V 



XIS 1 



3.8 

V 

Sw 5b 

00031280008 

pc 

2008-11-27 

1.9 

252 

V 



XIS 3 



3.7 

V 










HXD 


37.2 

2.6 

V 








Suz 6 

703022060 

XISO 

2008-12-02 

51.4 

1.7 

V 

Sw 6a 

00031280009 

pc 

2008-12-02 

0.3 

50 

X 



XIS 1 



1.8 

V 

Sw 6b 

00031280010 

pc 

2008-12-03 

1.7 

264 

V 



XIS 3 



1.7 

V 










HXD 


36.5 

2.3 

V 










\XMM 

\XMM 

1 ^ ^ ^ ^ ^ ^ ^ 

Suzaku/Swift (6,12) I 

Swift 1 


1 

1 RXTE 

1 1 

1 1 

1 1 

Ml 

2000 

2005 

2010 

2015 


Time [years] 

Fig. 1. Time-line including all observations of NGC 3227 discussed 
here. The blue region shows the time-range of sustained RXTE monitor¬ 
ing from 1999-Jan-02 to 2005-Dec-04 jUttley & McHardy|2005| l. The 
observations by Suzaku and Swift from the 2008 campaign are displayed 
in red, the times of the XMM-Newton and recent Swift observations are 
shown in blue. 


After creating calibrated event lists with filtered hot and bad pix¬ 
els, events in 10-12 keV are screened for enhanced rates due to 
particle flaring. As no pile-up is evident, we extract 0.3-10keV 
spectra from all counts within 40" of the central source position 
of NGC 3227. The same angular radius of the extraction region 
is used to extract the background from a position ~6'off-source 
on the same CCD-chip. 


2.3. Suzaku 


We extract both the data of the X-ray Imaging Spectrometer 
(XIS; Koyama et al.|2007| l and the Hard X-ray Detector (HXD; 
[Takahashi et al.|2007| l. We use data taken by the front (XIS 0,3) 
and back-illuminated (XIS 1) chips in the 3x3 and 5x5 edit¬ 
ing modes. We reprocessed the unfiltered event lists by applying 
the newest calibration available and screened the data with de¬ 
fault parameters. We then perform an attitude correction with 
aeattcor2 based on the 3x3 mode which comprises the bulk 
of data. The resulting spectra of the 3x3 and 5x5 modes are 
merged using mathpha for each XIS. Spectra were extracted from 
circular regions of ~93" radius and centered on the point source. 
We normalize all fluxes with respect to XISO. Our fits find that 
the flux normalization of the XISl and XIS3 spectra deviates 
by about 5% from that of the XISO spectrum, consistent with 
the Suzaku ABC guide (http://heasarc.gsfc.nasa.gov/ 


docs/suzaku/analysis/abc/), version 5.0. We perform si¬ 
multaneous fits of observations Suz 2-Suz 6 using the SimFit 
routines of ISIS (Kiihnel et al. 2015| l, and rebin their spectra 
to a combined minimum signal-to-noise ratio (S/N) of 25. Ob¬ 
servation Suz 1 is fitted separately and binned to a minimum 
S/N of 18. In both cases the binning does not exceed the ener^ 
resolution of the XIS detectors of 150 eV at ~6keV in 200® 
Due to calibration uncertainties caused by the Si- and Au-edges 
we exclude data in the 1.72-1.88 keV and 2.19-2.37 keV energy 
bands. For Suz 1 the first interval is extended to 1.5-1.88 keV 
because of further insufficient calibration just above 1.5 keV. 

The non-imaging HXD-PIN data are extracted for the whole 
field of view of 34' x 34'. For a list of total counts and good 
exposure times see Table The flux normalization of the HXD 
with respect to XISO has a large uncertainty because of the low 
signal to noise ratio of the HXD data, but is consistent with what 
is expected for the current calibration. We therefore fix the HXD 
flux normalization constant to its nominal value of 1.16. 


2.4. Swift 


The X-ray Telescope (XRT) on board Swift ( [Burrows et al. 2005 1 
provides imaging and spectroscopic capabilities for the energy 
range 0.5-10 keV. Except for observations Sw 2a,b and Sw 3a, 
which were taken in windowed-timing (WT) mode, all observa¬ 
tions were performed in photon-counting (PC) mode. For the PC 
mode we allow grades from 0 to 12. Only grade 0 events are ex¬ 
tracted for the WT mode. We extract source counts within 30' of 
the source position. The background is extracted from circular 
regions of at least 60" radius on a position on the chip which is 
free of background sources. We exclude observation Sw lb from 
our analysis because of an excess above 7keV of unknown ori¬ 
gin. Observation Sw 6a was also ignored due to an insufficient 
number of detected counts. 

The Ultraviolet and Optical Telescope (UVOT) on board 
Swift observes with up to six filters (Roming et al. 2005 UVW2, 


^ https://heasarc.gsfc.nasa.gov/docs/astroe/prop_ 
tools/suzaku_td/nodeIQ.html 
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UVM2, UVWl, U, B, and V). The UVOT image data have been 
summed using uvotimsum (V. 24Jan2014_V6.15.1). A circular 
region of 5" around the source position was used for the extrac¬ 
tion of counts of the central core component. The background 
was defined as an annulus around the source with 13" and 26" 
inner and outer radius, to minimize contributions from the host 
galaxy. The resulting files were converted into pha files using the 
uvot2pha task and de-reddened as discussed in Sect. 13.41 


3. X-Ray Spectral Analysis 

The main focus of the data analysis lies on the simultaneous 
Suzaku and Swift data describing an additional absorption event. 
In the following, we first build up a spectral baseline model us¬ 
ing the Suzaku data. With help of this model, we re-visit the two 
archival XMM-Newton observations and finally model each set 
of Suzaku and Swift observations simultaneously. 


3.1. A baseline model based on Suzaku data 

Before we start to develop a baseline spectral model for 
NGC 3227, we briefly review the previous attempts to de¬ 
scribe the spectrum. Lamer et al. ( 2003| l describe the 2000 
XMM-Newton observation with a power law absorbed by neutral 
or lowly-ionized material (hereafter WAj) and an unabsorbed 
power law. The slope of the absorbed power law was difficult to 
constrain. As shown by Lamer et ar] ( |2003] l, this model is highly 
degenerate in f (two orders of magnitude) and Ah. We And that a 
lower degree of degeneracy can be obtained by assuming a par¬ 
tial covering scenario which yields comparable statistics. The 
lack of sufficient RGS data and the relatively low signal-to-noise 
did not allow [Lamer et al.|P003| l to constrain more absorbers. 

In contrast, the 2006 XMM-Newton observation caught the 
source in a typical hard X-ray flux and spectral state, i.e., unab¬ 
sorbed by moderately Compton-thick gas. The high soft X-ray 
flux and the long exposure time allowed [Markowitz et al.| ( |2009| l 
to obtain a high S/N RGS spectrum. They constrained two layers 
of ionized absorption, one intermediately ionized, log^ ~ 1.45 
and one highly ionized, log^ ~ 2.93 absorber. We call these 
WA 2 and WA 3 in the following. [Markowitz et al.| ( |2009| l also de¬ 
tected an absorber with a small neutral column, which is similar 
to WAi seen in 2000 by [Lamer et al.[ ( |7003[ |. 

In deriving the baseline model for our analysis, we use the 
high signal to noise Suzaku data and assume from now on that 
WAi, WA 2 , and WA 3 - assumed to represent physically distinct 
absorbing media with different levels of ionization - are present 
in all XMM-Newton, Suzaku, and Swift observations. The as¬ 
sumption that WAI, WA2 and WA3 are present in all observa¬ 
tions and stable over durations of many years is a simplifying 
one; in some cases, a given absorber may not be statistically 
required in the fits, particularly if the column densities of the 
other absorbers are such that they dominate the spectral shape. 
In addition, as demonstrated below, we observe a wide range in 
column densities in WAI and WA2 across spectra taken several 
years apart. A simple assumption is that each absorber remains 
present in all observations, and has varied primarily in column 
density, and not ionization state. However, we cannot exclude the 
possibility of some observations capturing an additional physi¬ 
cally distinct absorber with similar ionization level moving into 
the line of sight. 

Figure [^ shows the steps followed to And the best-fit base¬ 
line model on the 2008 Suzaku spectra motivated by Lamer et al. 
( [2003 | l and [Markowitz et ar] ( [2009l l. We use observation Suz 3 
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in this example, which has the best combination of count rate 
and strong absorbing column density. The of each fit is la¬ 
beled in the figure and shows the gradual improvement of the 
model. In Fig. [^ the high-energy X-ray power law (HXPL) ab¬ 
sorbed by one layer of neutral absorption (CA). Apart from the 
prominent Fe Ka line around 6.4 keV, the strong residuals below 
4 keV suggest a partial covering scenario similar to that seen in 
the 2000 XMM-Newton data. Assuming partial absorption by a 
neutral medium results in the residuals shown in Fig. Be¬ 
low 3 keV the residuals still show strong positive and negative 
excesses. Replacing the neutral absorber by a warm absorbing 
layer of very low ionization results in substantially flattened 
residuals (Fig.j^). We identify this component with WAi. 

While this model is able to describe the smooth turnover of 
the partial coverer, some residuals remain. Based on the 2006 
XMM-Newton observation, we model these residuals by includ¬ 
ing two additional layers of ionized absorption, WA 2 and WA 3 
(Fig. |2j l). Here, following the RGS analysis of |Markowitz et ah] 
POOi^ the absolute redshift of WA 3 is fixed at -0.00302. The At 
requires the intermediately ionized absorber WA 2 (log^ ~ 0.5) 
to be dominant and partially covering. 

Leftover residuals in the final continuum fit are line-like fea¬ 
tures in the soft-band around 0.9 keV, which are most likely due 
to Neix emission, and the strong signature of a Fe Ka/Fe Kf 
complex between 6.4 and 7.1 keV. These soft and hard compo¬ 
nents can either be modeled individually using empirical emis¬ 
sion lines, or by adding ionized reflection (modeled with xillver; 
[Garcia et al.[[2013} , where we require that the reflected contin¬ 
uum has the same shape as the primary continuum. Both ap¬ 
proaches describe all of the excess components equally well and 
without appreciably changing the continuum parameters. Be¬ 
cause of the unique physical interpretation and the smaller num¬ 
ber of free parameters, we continue the modeling by using xil¬ 
lver. The resulting model describes the data very well (Fig.j^). 

3.2. Re-Visiting the XMM-Newton spectra with the baseline 
model 

We now test whether the baseline model also describes the 
XMM-Newton observations that have been used to motivate this 
model and its absorbers. The spectra and model components 
of the 2000 and 2006 XMM-Newton observations are shown in 
Fig . The continuum is described by an absorbed high-energy 
power law, leaked emission, and unblurred ionized reflection 
components. Table [^and Table [^list the best-fit parameters. 

Because of the inability to constrain the covering fraction 
of the mildly (WAi) and highly (WA 3 ) ionized absorbers, we 
always require them to fully cover the compact source, while 
their column densities and ionization parameters are all left free 
to vary. 

The continuum as well as the Fe K ajp lines of the 2006 
XMM-Newton observation are well fitted with a combination of 
an incident power law and ionized reflection component ( [Garcia 
[et al. 2013| xillver). The two modest absorption dips at around 
0.75 keV and 0.9 keV require WA 2 to have log^ ~ 2 and to cover 
the X-ray source by ~50%. Figure [^ shows the 0.5-2.3keV un¬ 
folded spectra of both the Suz 1 observation from 2008 and the 
XMM-Newton observation from 2006. The two absorption fea¬ 
tures enclosed by the gray shaded region are found at about the 
same energy in both spectra, arguing for WA 2 as a common ori¬ 
gin. The similar spectral shape prefers WA 2 to be partial cover¬ 
ing in both observations. In contrast to the 2000 data, a steep soft 
excess is seen in 2006, which we describe by a soft power law 
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Fig. 2. Example of the steps followed to find the best-fit baseline model 
on the Suzaku spectra (we use observation Suz 3 in this example), a 
Fit to the spectrum with a single power law absorbed by neutral gas. b 
The neutral absorber is partially covering the incident power law. c Due 
to characteristic residuals we replace the solely neutral absorber by an 
ionized one and d and two more warm absorbers leading to a stratified 
composite of differently ionized absorbers, e The remaining residuals 
in the soft band and around 6.4 keV are flattened when finally adding a 
reflected power law in panel. 


(SXPL), consistent with the results of [Markowitz et aLl ( |2009| l. 
The fit converges to a statistic of;if^/dof = 1851/1542. 

While the data from 2006 are described well with the partial 
covering absorber WA 2 , the 2000 data require a dominant partial 
covering column of lowly-ionized gas (WAi). Given the strong 
X-ray absorption, it is not clear if a soft excess is present in the 
spectrum. If our partial-covering model is correct and the contin¬ 
uum below ~2 keV is dominated by “leaked” emission, then the 
soft excess (T ~ 3) must have had a negligible presence. Given 
the low signal-to-noise ratio of the data, with the exception of 
the column of WA 2 the parameters of the two highest ionized ab¬ 
sorbers WA 2 and WA 3 cannot be constrained, and are assumed 
to be equal to those found in the 2006 data. The emission fea¬ 
ture at ~0.88keV is described by a narrow unresolved Gaussian 
for Fe Lxvii at 0.826 keV and an O viii radiative recombination 
continuum (RCC) at 0.871 keV. We find the best-fit statistic to 
be;^'^/dof = 901 /847, which is consistent with the one found by 
[Lamer et al.| ( |2003| ). 

We conclude that the XMM-Newton data of both the ab¬ 
sorbed ( 2000 ) and relatively unabsorbed (2006) observations of 
NGC 3227 can be described with our baseline model regard¬ 
ing the absorber structure and unblurred reflection. The spectral 
variability between both is dominated by the mildly ionized ab¬ 
sorption component WAi. 

3.3. Simultaneous fit to all Suzaku spectra 

The success of the baseline model in describing both the individ¬ 
ual Suzaku data and the XMM-Newton-dsLtn suggests that we can 
use it to model all 2008 Suzaku spectra simultaneously. The fit 



Fig. 3. The 2006 (top) and 2000 (bottom) XMM-Newton spectra, the 
best fit model components, and the complete model shown with thick 
and thin solid lines, respectively. Both spectra are described by the base¬ 
line model that comprises a high energy X-ray power law with a reflec¬ 
tion component and absorption by three layers of lowly, intermediately 
and highly ionized gas as well as a Galactic column. The source is par¬ 
tially covered by WA 2 for the 2006 and by WAi for the 2000 obser¬ 
vation. The data from 2006 also require a steep soft excess. The sharp 
residuals at '-1.85 keV are likely due to calibration uncertainties around 
the Si K edge. 



2 _I_^_I_I_I_^^_I_^_I_^^_ flhffiMil 1i 

0.5 1 1.5 2 

Energy [keV] 


Fig. 4. Comparison of the 2006 XMM-Newton and Suz 1 observations 
in the 0.5-2.3 keV range. The gray area marks the energy range com¬ 
prising two modest absorption features that are due to an absorber of 
similar i.e., WA 2 . 
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Table 2. Continuum parameters of the spectral fit of the baseline model 
to the relatively unabsorbed 2006 XMM-Newlon observation. Fixed pa¬ 
rameters are marked with an asterisk. We assume that the incident hard 
X-ray power law (HXPL) of the ionized reflection and the continuum 
power law are identical. To reduce model degeneracies, we freeze the 
photon index to the value found by [Markowitz et JI]p009[ l, who find 
WA 3 to be outflowing. We therefore also freeze its absolute redshift to 
the according value. The soft excess is modeled by a power law (SXPL). 
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results indicate that the parameters of WAi and WA 3 stay con¬ 
stant during the 2008 observational campaign, while the spectral 
variability is dominated by WA2. We have to model Suz 1 sepa¬ 
rately, as it does not share some of the time-independent param¬ 
eters. Table|4]lists the results of the fit. 

In more detail, the time-independent parameters include the 
detector constants, the iron abundance Zpe of the ionized reflec¬ 
tor (a simultaneous fit of Suz 2 to 6 yields Zpe = 2.81 +0.17) and 
all parameters of the non-varying ionized absorbers WAi and 
WA 3 . The mildly and highly ionized absorbers WAj and WA 3 are 
assumed to fully cover the central source. WAi has a relatively 
low column density of 0.137^QgQg x 10 ^^ cm“^ of mildly ionized 
(log^ = -0.29+°;°^) gas for Suz 1 and 0.068+^25 x lO^^cuT^ 
(log^ = -0.9 + 0.6) for the later 5 observations. While both 
log ^ are consistent with each other, the column found in Suz 1 
is slightly higher. For WA 3 the column densities are consistently 
around ~ 4 x lO^^cm^^, while the ionization parameter differs 
between the single fit of Suz 1 (log^ = 3.44;!;^^) and the other 
observations, which have (log^ = 4.17;|;g22)' Figure shows 
how both instances of WA 3 affect the high energy power law. 
The model with log^ ~ 3.4 for Suz 1 imprints a clear absorption 
feature around 0.9-1.0keV due to Neix and Nex. A column of 
comparable ionization was also constrained by [Markowitz et ah] 
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Table 3. Continuum parameters of the spectral fit of the derived baseline 
model to the absorbed 2000 XMM-Newton observation. Fixed parame¬ 
ters are marked with an asterisk, parameters adopted from the fit to the 
2006 XMM-Newton observation with f. The photon index of the reflec¬ 
tion component is tied to the one of the high energy power law. The 
model also includes an unresolved Gaussian line at Ec = 0.826 keV 
identified as Fe L xvii and a O vin radiative recombination continuum 
(RRC) with an edge energy of 0.871 keV. 
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Fig. 5. Effect of WA 3 if it were to fully cover the background power law. 
Shown in red is the warm absorber of observation Suz 1, black shows 
the effect of the absorber found in the remaining Suzaku observations. 
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MJD (start of observation) 


Fig. 6. Temporal evolution of the parameters describing the variable 
intermediately ionized absorber WA 2 and the luminosity in the 0.6- 
lOkeV X-ray band for the six Suzaku ohservations. 



12 5 10 

Energy [keV] 


Fig. 7. Effect of the partial coverer WA 2 onto a power law (thick lines 
with line styles given in the figure) as part of the best-fit model (Eig.[^. 
The thin and solid lines of lighter color show the effect the coverer 
would have if it were to fully cover the source. 


( |2009| l but disappears in Suz 2 to Suz 6 . It may, however, still 
be intrinsically present but undetected because of degeneracies 
within the model description. The highly ionized absorber WA 3 
(logf ~ 4) shows weak absorption lines from highly ionized 
species of ion (Fig.[^. These are difficult to constrain in a single 
Suzaku observation but are consistent with the detailed study of 
Suz 1-Suz 3 by |Gofford et al.| ( [2013| ), who show the existence of 
absorption lines due to He-like and H-like Fe. 

The time-dependent group parameters of the model repre¬ 
sent the effects of variability. In particular we mention the vari¬ 
ability of the intermediately ionized partially covering absorber 
WA 2 . Simultaneous fitting isolates WA 2 as the only variable 
layer of absorption within the complex interaction of the three 
absorbers contained in the model. Figure shows the time evo¬ 
lution of the fit parameters Ah, log /cvr and the X-ray luminos¬ 
ity Lx. The covering fraction scatters around ~0.9. Only Suz 1 
shows a lower value of ~0.74. The column density peaks in a 
complex fashion between 5-18x10^^ cm^^ across all six obser¬ 
vations. The ionization parameter first decreases from log^ ~ 2 
to nearly zero and then increases again to values of log f ~ 1 
for the last two observations. In Fig. |7]we illustrate the effect of 
WA 2 onto the power-law continuum if this absorber had a cover¬ 
ing fraction of unity (see also Fig.[^. This representation of the 
model helps to identify the most prominent contributing absorp¬ 
tion features below 1 keV, which are mainly due the He-like and 
H-like ions O vii, O vm, Ne ix, and Ne 

Other variable parameters are the photon-index F that scat¬ 
ters between ~1.6 and ~1.7. We also find intermediate ioniza¬ 
tion states of the reflecting and absorbing material ranging from 
log^ ~ 0.8 up to ~2. 

The results of the fits are shown in Figure We show the 
XIS and HXD spectra of all observations including the best- 
fit models. The residual panels suggest a fit that describes the 


^ Note that because in reality the warm absorber is only partly covering 
the source, the fit is not particularly sensitive to these narrow features 
and mainly driven by the characteristic curvature below 2 keV caused 
by the intermediate ionization of WA 2 . 
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Fig. 8. Spectra and residuals in units of standard deviations for a fit of 
the baseline model to the Suzaku spectra of the 2008 campaign (Ta¬ 
ble]^. A single fit is performed for the high state observation (Suz 1) 
while the remaining observations (Suz 2-Suz 6) are fitted as part of one 
simultaneous fit with tied time-independent parameters. The statistics 
printed in the residuals panels correspond to the best-fit on each obser¬ 
vation, using the appropriate time-dependent parameters. 


continuum well with adequate statistics. Figure in the ap¬ 
pendix gives an overview of all Suzaku spectra shown together 
with the full model and the best-fit model components. Clearly 
visible is an emission excess above the power law in the en¬ 
ergy range covered by ^Mzaku/HXD/PIN, most likely due to the 
Compton hump peaking around 30keV. The 6.4 keV Fe Ka line 
is self-consistently modeled by reflection off Compton-thick gas 
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(xillver), as is the extra emission below IkeV. Figure [^il¬ 
lustrates that the reflection parameters are mainly derived from 
that soft emission and the iron line complex as dominant fea¬ 
tures. The HXD/PIN data quality is good enough to constrain the 
normalization of xillver. No extra soft emission component is 
needed in contrast to the 2006 XMM-Newton observation. Al¬ 
though we can not rule out its presence, the soft X-rays are well 
modeled by a combination of ionized reflection and leaked coro¬ 
nal emission described by the HXPL. The (narrow) Fe Kcr line 
has an equivalent width of EW ~ 130eV and is well described 
by xillver. No further line emission is required, indicating that 
the line is consistent with being completely due to reflection off 
Compton-thick material. 

The complexity of our model leads us to study possible de¬ 
generacies by calculating contours for all pairs of time- 
dependent parameters. These contours show that the simultane¬ 
ous fit is remarkably robust. As an example. Fig. [^ shows the 
contours for observation Suz 2. The overall correlations are sim¬ 
ilar for the other observations but differ in extent depending on 
the data S/N. A mild degeneracy between the covering factor 
with both the normalization and the photon-index of the high 
energy power law is present, as are modest degeneracies be¬ 
tween the ionization state and the column density of the WA 2 
as well as the photon-index and the normalization of the power 
law. The contours of xillver show an overall degeneracy with 
most other parameters, which can be removed by assuming that 
photon-index is a time-independent parameter. Finally, turning 
to the time variability of the absorber. Fig. shows the con¬ 
tours of all time-dependent parameters related to Ah for all ob¬ 
servations and in particular reveals a clear variability in column 
density across the observations. Note that the contours includ¬ 
ing the normalization and ionization of the unblurred reflection 
component xillver reveal no significant variability. 

3.4. Analysis of Swift XRT data 

For each Suzaku observation the archive contains two simultane¬ 
ous Swift observations. 

We model 2008 Swift XRT simultaneously with the same 
baseline spectral model but independently of the simultaneous 
fit of the Suzaku data. Tests show that the Swift data are consis¬ 
tent with the Suzaku data and do not contribute with additional 
information due to the lower effective area. We therefore impose 
constraints gained with Suzaku for certain parameters that can 
neither be constrained with Swift nor identified as variable with 
Suzaku. All derived parameters are consistent with those found 
for Suzaku data. 

We also analyze the most recent Swift observations from 
2013-2015. The data are well described with a single power 
law absorbed by both Galactic and intrinsic, neutral absorption 
modeled with tbnew. We find a very low degree of absorp¬ 
tion (Ah < 0.34 X lO^^cm^*) consistent with the predominant 
relatively unabsorbed state during the RXTE monitoring over 
6.9 years. 


4. Multiwavelength Data Analysis 

4 . 1 . NIR-X-ray SED and reddening 

The broadband SED (Fig.[TT]i is based on Swift/UWOT data from 
the observation Sw 4a and allows us to estimate the amount of 
reddening due to dust in the line of sight and to draw conclu¬ 
sions on the dust content of the X-ray absorbing cloud. Model 
degeneracies, however, require us to also include data from other 
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Fig. 9. Color coded contours for relevant spectral parameters of 
Suz 2. The 68.27%, 90% and 99% contours correspond to the solid, 
dashed and dotted lines. The color code ranges from brown (small Ax^) 
up to dark blue (large Ax^). 


wavebands. Those are photometric data from ground- and space- 
based facilities of different spatial resolutions and apertures: 
WISE ( [Wright et al.l20T0| and 2MASS ( [Skrutskie et ah] 


20061 


data from the All-Sky Source Catalogs, data ( [Ramos Almeida| 
[et al.||20TT) from the NASA 3 m IRTE telescope (NSFCam at 
L-band), the 3.8 m UKIRT telescope (IRCAM3 at M-band), the 
Gemini-South telescope (T-ReCs, N-band) as well as HST NIC- 
MOS data in the Alters E160W and E222M ([Quillen et al.|1999|l, 
HST EGS data at E550W (jBentz et al.|2009t and HST ACS data 


at 3300A ( Munoz Marin et ^.[2007 1 . This last data point, taken 


in 2002, is roughly consistent with the UVOT flux taking into 
account the different spatial resolution as well as likely long¬ 


term source variabilitjy In addition, we simulate Spitzer data 
based on the AGN contribution to the IR, as modeled by |Hernat>] 
[Caballero et al.|p015] l. A cross-calibration factor between WISE 
and Spitzer of 1.6 is consistent with flux calibration uncertainties 
and aperture effects ( [Hernan-Caballero et al.|20T5) . Eor model¬ 
ing the continuum, we only take into account data with PSEs that 
are comparable to the one of 5w//f/UVOT, i.e., WISE, 2MASS 
data with PSEs of 5"-10"radius, as well as SwiftPiKT data. Sys¬ 
tematic effects and model degeneracies due to aperture effects of 
the instruments and a lack of simultaneous data lead us to ab¬ 
stain from performing a statistical minimization at energies 
lower than the X-rays. Instead, we use sub-arcsec Spitzer, HST 
and the ground-based data that are only sensitive to the AGN 
core to constrain the AGN-related model components, i.e., the 
dusty torus and the accretion disk emission in the IR and UV, 
respectively. While the Swift data provide the only variability 
information on timescales of weeks, the X-ray band turns out 
to be the only band with strong variability and we can use the 
already known model between 0.6-10 keV. The IR-UV bands, 
in contrast, show only minor variability up to 25% on the time- 


The extended flag in the ALLWISE catalog introduces problems with 
the profile-fit photometry, which is why we add 10% systematic uncer¬ 
tainties to 2MASS and WISE data as well as 5% for UVOT, which has 
a similar PSF. 
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Fig. 10. Contours relating all time-dependent parameters of all observations to Ah of WA 2 . The solid, dashed and dotted lines represent the 68%, 
90%, and 99% confidence contours. Observations Suz 1 to Suz 6 are both marked with numbers and color-coded in black, blue, green, red, orange, 
and purple, respectively. 
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IQia ;^q14 ^q15 ;^q 16 jgl8 
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Fig. 11. Example-SED of NGC 3227 illustrating the coverage from 
the MIR (WISE: red circles; Spitzer. purple circles and HST\ orange 
triangles) over the NIR (ground-based data: brown squares), the opti¬ 
cal (2MASS: blue circles; HST) and UV (Swi/r/UVOT: turquoise cir¬ 
cles) to the X-ray band up to lOkeV (Swi/I/XRT: green circles). The 
UV and X-ray data were simultaneously measured at 2008-11-21 while 
the 2MASS and WISE data are archival catalog data. The Spitzer data 
are simulated data based on the AGN contribution (|Heman-Caballero| 
[eTaLjllOTS) . Other data are non-simultaneous photometric measure¬ 
ments from the literature. See text for details. The dashed lines corre¬ 
spond to the unabsorbed and unreddened model components, the thick 
solid line to the overall reddened and absorbed model. 


scales monitored by SwiftfUWOT. In the following we describe 
our best fitting model for data below the X-rays, for which we 
find;rVdof = 125/23. 

The reprocessed MIR/NIR emission probed by Spitzer, 
ground-based and HST data at 1550A and 2300A have been 
fitted previously with clumpy torus models by, e.g., [Ramos | 
[Almeida et al.| ( |2()ll| l or [Alonso-Herrero et al.[p011[ l. We in¬ 
stead use a phenomenological description of three black-bodies 
from 145 K up to 1480 K that is supported byjCalderone et al.[ 
( | 2 ^ . 


We take account of the host galaxy starlight using a SWIRE 
template for a Sa spiral galaxy. The broad-band template is 
co nstructed to match a stellar template spectrum (3500-7000A) 


byfWinge et al. 

1995 

1 at 5500A for a normalization of one. Be- 

easures extra emission that can not be fit- 
disk. According to [Davies et al.[ ([2006|l, 

nilano.inaf.it/~polletta/templates/ 

low 4000A, UVOT m 
ted with an accretion 

^ http: //WWW. iasf-r 

swire_templates.html 


the central 30 parsecs show remainders of starburst emission ac¬ 
counting for 20% to 60% of the galaxy’s bolometric luminosity. 
We therefore add a starburst template ( Kinney et al.|1996| l below 
4000A. This is also confirmed by an independent photometric 
measurement of the starlight by HST at 3300A (lower triangle 
in Fig. 11 1 that perfectly matches the combined template for a 
normalization of one. We find a template normalization of ~4, 
which is due to uncertainties in the derivation of the stellar tem¬ 
plate and aperture effects. 

The accretion disk emission is described with a multi¬ 
temperature disk-blackbody (diskpn, ~ 1 x 10^ K at 

fin - 6rg reddened with a column of 1.2 x 10^'cm“^ or 
Ay ~ 0.45 mag). Its normalization is constrained by HST data at 
5470A and 3300A. The data have been taken 2 years apart from 
each other, which can explain their scatter around the model. 
[Koratkar & Blaes[ ( |1999| l confirm that the value found for r^ax is 
consistent with sample studies of quasars with minimal extinc¬ 
tion as well as the expected inner disk temperature of a thin a- 
disk. Degeneracies between the applied reddening and the inner 
disk temperature are reduced due to the consistency of the red¬ 


dening column with results from the literature. According to Ko¬ 
ratkar & Blaes ( 1999[), one usually assumes extinctions of E(B- 


V)~0.05-l mag, (Ah ~ 0.04x 10^^-0.8x 10^^cm ^foraGalac- 


tic gas-to-dust ratio). According to Komossa & Finkj ( [1997| l; 
[Kraemer et al.[ ( [2000| l; [Crenshaw et al.[ ( [2001| l, NGC 3227 has 
shown a rather large reddening column of ~ 0.2 x 10^^ cm“^ that 
is claimed to be consistent with a persistent, dusty, “lukewarm” 
absorber at ~ lOOpc distance. Swift observations from 2013- 
2015 likely measure the same distant absorber but with only half 
the column, which is consistent with the reddening found in our 
model and maybe also by Winge et al. ( |1995| l. 

Compared to the variable X-ray-absorbing column densities 
measured with Suzaku and Swift, the amount of UV extinction 
is associated with a column that is ~ 70 times smaller. Even for 
larger reddening and therefore larger inner disk temperatures, 
the severe mismatch to the X-ray column persists. This is unex¬ 
pected if both the X-ray and UV-absorber would be of the same 
origin. 


For the extinction at optical wavelengths (Ay; Fitzpatrick 

om Crensha' 
int according to 


1999 1 we use the reddening-curve c(/l) from Crenshaw et al. 
( |200i[ l that attenuates the intrinsic flux F\ 

^’obs = ^’int ■ (2) 

with the extinction coefficient at the wavelength A, = 
c(/I)£’(B - V) = c{A)Ay/Ry. We choose Ry = 3.1 and estimate 
the equivalent column density via the assu med Galactic gas-to 


dust ratio Ah = Ay -2.69x10"^* cm ^mag ^ ( [Nowak et al.|2012 1. 


Note that the standard Galactic reddening curve does not have to 
apply for the intrinsic absorber in NGC 3227 where dust grains 
might have a different composition and size distribution com¬ 
pared to the Galaxy ([Crenshaw et al.|2001[). 
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Fig. 12. a Light-curve of the unahsorbed X-ray luminosity in the 0.6- 
lOkeV band; b light-curve of monochromatic luminosities in the B- 
band (blue circles), R-band (red triangles) and I-band (brown squares) 
for the 84 cm and 1.5 m SPM telescopes; the R-band and I-band lu¬ 
minosities are normalized to the weighted mean of the B-band with 
relative factors of 0.72 and 0.32, respectively; c/d degree of opti¬ 
cal polarization p^pt (black squares) and electric vector position angle 
(PAopt, blue triangles) for the B/R-filters. The I-filter did not provide 
polarimetry data. The central wavelength of the B/R/I-filters are given 
as 0.438/0.641/0.798 pm by |Bessell et aLl ( [T9^ . 


4.2. Optical Polarimetry 


A portion of the optical AGN emission undergoes scattering into 
our line of sight resulting in a few percent of linear polariza¬ 
tion. Investigations prefer the scenario of polar-scattering for 
NGC 3227 ( [Smith et al.||2008] l. This result implies an interme¬ 
diate inclination of the nucleus towards the line of sight, simi¬ 
lar to what is predicted by clumpy torus models that are tested 
for NGC 3227 in Sect. |5.5| We therefore investigate additional 
optical polarimetric data coincident to the 2008 SwiftjSuzaku 
campaign that have been taken between 2008 October 26 and 
2008 December 03 with the 84 cm and 1.5 m telescopes of the 
San Pedro Martir observatory (SPM) in Mexico. The data were 
taken with the B hlter (Johnson system) as well as R and I biters 
(Cousins system) according to the standard photometric system 
( Bessell| 20051. For the magnitude-bux conversion we use the 
zero-point magnitudes from|Bessell et al.|(|1998|l. 


We compare the optical photometric luminosities with the X- 
ray luminosities in Fig. [T^ The relatively sparse sampling pat¬ 
tern precludes a detailed cross-correlation analysis between the 
X-ray and optical light curves. The two bottom panels show po¬ 
larimetry data for the B- and R-band. The degree of polarization 
p and the electric vector position angle (PA) scatter around a 
few percent and 100°-150°, respectively, which is roughly con¬ 
sistent with the results from [Smith et al.| ( |2008| l. They bnd that 
the PA is constant over wavelength with a value of PA ~ 135°. 
In contrast we observe different weighted averages between the 
B-band (PA ~ 140°) and R-band (PA ~ 119°). [Mundell et al.[ 
( [1995a[ l bnd both signatures for a narrow line region in high- 
excitation O III emission lines (PA ~ 30°) and two compact radio 
cores along a PA of ~170°. [Mundell et al.[ ( [r995b[ ) identify a 
much more extended radio source with PA ~ 158° that shows 
a rotating disk-like structure on scales of a typical host galaxy 
with an inclination of ~56°; any connection to the AGN on much 
smaller scales is not immediately obvious, especially given the 
uncertainty regarding the AGN system inclination. We support 
the suggestion of [Mundell et al.[ ( |l995a[ ) and [Smith et al.[p008[ l 
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Fig. 13. Evolution of the column density Ah of WA 2 over time as found 
in the Suzaku (black triangles) and Swift (gray circles) data. We plot 
uncertainties according to the 90% conbdence level. The two data points 
derived by Swift that are cut off by the plotting window have the values 
Ah = 42/32 ^ lO^^cm^^ and Sd/jg x lO^^cm^^. 
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that the cone of excited O in gas may be oriented along the sym¬ 
metry axis of the AGN. The angle enclosed between the cone 
and the measured average PA in the B-band is therefore ~102°, 
which is close to a PA of 90°. Together with the degree of polar¬ 
ization we therefore bnd the optical emission in NGC 3227 to be 
consistent with polar-scattering, similar to what is suggested by 
[Smith et al.[ ( |20^ . 

When facing additional systematic uncertainties in the po¬ 
larimetry data we can not conclude any signibcant variability. 
The optical luminosities, in contrast, show less than 40% vari¬ 
ability with systematics most likely affecting the I-band data 
point at MJD 54796.7. 


5. The properties and origin of the absorbing gas 

5.1. Summary of observational results 


The 6.9 years of sustained RXTE monitoring from 1999 Jan to 
2005 Dec ( [Rivers et al.[[20lTb i caught two discrete eclipses. 
The brst ~80 d long event in 2000 was analyzed by[Lamer et al. 


P003[ l and is dominated by a lowly-ionized absorption compo¬ 
nent (log^ ~ 0.4; WAi in our model). The second event lasted 
2-7 days and was also lowly- or at most moderately-ionized 
( [Markowitz et al.|2014[ l. Based on this monitoring, the relatively 
unobscured observation by XMM-Newton in 2006,, and recent 
Swift observations from 2013-2015, NGC 3227 was usually rel¬ 
atively unobscured before 2006 and after 2013. Thanks to these 
previous observations and in particular the 2008 Suzaku/Swift 
campaign, we can build a clearer picture of the behavior of 
the variable absorption components in NGC 3227 over the last 
15 years and also address the properties and origin of the ab¬ 
sorber during the 2008 campaign. 

During the 2008 campaign, the observed spectral variabilit 
was dominated by absorption due to intermediately-ionized ga: 
(log^ ~ 1; WA 2 in our baseline model) with a time-variable 
complex and irregular density prohle, and columns spanning 
~5-16 X lO^^cm^^ (Fig. 131. The column densities derived by 
Swift are mostly consistent with those from Suzaku within 2cr. 
The event duration is > 35 days. Due to the lack of concurrent 


® We caution that WA 2 as modeled during the 2008 campaign and WA 2 
as measured during the 2006 XMM observation may represent different 
components: the ionization parameters may be similar but the column 
densities differ by more than a factor of 50. 




























































T. Beuchert et al.: A Variable-Density Absorption Event in NGC 3227 mapped with Suzaku and Swift 


100 

80 
u 

I 60 

M 
u 

QJ 

^ 40 

20 

0 

2x10'*^ 3x10^^ 4x10^^ 

Lx/Nu [ergs-^cm^] 

Fig. 14. Correlation between the ionization parameter ^ and the X-ray 
luminosity Lx as a proxy for the ionizing luminosity. The numbers give 
the Pearson and Spearman rank correlation coefficients as well as the 
null hypothesis probabilities. 
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RXTE monitoring in 2008, we have no firm information on cloud 
ingress/egress. 

5.2. Luminosity and ionization 

We now proceed to study further the nature and properties of 
WA 2 by testing if there is any relation between its ionization 
state and the measured source properties. 

Because of Eq. 1, we expect a linear relation between f and 
the ionizing luminosity. We first assume that the measured X- 
ray continuum is a direct proxy for the variation of the ionizing 
continuum, i.e., we assume that the measured X-ray luminosity 
by Suzaku, taken here from 0.6 keV to 10 keV, is directly propor¬ 
tional to the overall ionizing luminosity Lion- Figure shows 
that there may be a trend in the data that higher luminosities 
also imply larger The Pearson correlation coefficient for f as 
a function of Lx is rp_x = 0.70 with a null hypothesis probability 
of Pp,x - 0.12. However, the Pearson correlation coefficient is 
sensitive to extreme outliers, i.e., the data from Suz 4. As a sec¬ 
ond check we also calculate the Spearman rank correlation coef¬ 
ficient, which is less sensitive to extreme outlying points than the 
Pearson coefficient. We find rs,x = 0.71 {Ps,x - 0.11). The data 
thus do not allow us to claim a direct relation between ^ and Lx', 
any correlation is tentative at best. More observations, especially 
at lower luminosity states, are needed. 

5.3. The iocation of the absorber 

Here we provide estimates on the distance of the variable ab¬ 
sorber detected in the 2008 campaign from the SMBH. Through¬ 
out this section we make the simplified assumption that the 
mapped column density profile is due to a single spherical cloud 
of uniform density, illuminated by a central source and there¬ 
fore ionized. For the ionization state and column density we will 
use the average fit parameters from the Suzaku campaign, i.e., 
log^ ~ 1.1 and Ah ~ 11.2 x lO^^cm^^. 

Many of the properties of the cloud also depend on the ion¬ 
izing luminosity. Lion- As we showed above, there is only a ten¬ 
tative relation between the ionization parameter and the X-ray 
luminosity as proxy for the ionizing luminosity, i.e., the ioniza¬ 
tion state of the absorber is mainly determined by the source be¬ 


havior in the UV, which is inaccessible to our measurements and 
hidden in model degeneracies. We can place a constraint on the 


source’s UV emission from the measured SED shown in Fig. 11 


In order to derive the ionizing luminosity between 13.6eV 
and 13.6 keV (1-1000 Ryd), we assume that the UVOT data 
are non-variable within the 2008 campaign. We find Lion = 
8.9 X lO'*^ erg s ' for the sum of the non-variable disk-blackbody 
and the average of the X-ray power laws of all observations. 
The average bolometrrc luminosity is measured as Lboi = 1.3 x 
lO"*^ erg s“*. Using the assumed black hole mass, we find an Ed¬ 
dington ratio of TEdd ~ 0.005. 

The measured value of Lion is highly model dependent. 
We therefore compare our results with those of |Vasudevan| 
& Fabian ( 2009|l and Vasudevan et al. p010|l, who present 


strong sample studies for estimating the energy output of AGN 
but lack extended data coverage. They show that the ratio of 
the UV-luminosity below lOOeV to the bolometric luminos¬ 
ity, i'0°0i36-0.ikev/^Boi = 0.21... 0.59 for values of the Edding¬ 
ton ratio /ledd = pBoxIPpid ranging from 0.01 to 0.61. In con¬ 
trast to our study, |Vasudevan et al. ( 2010| l include nuclear IR 


that r'™ 

Liiai i-o.0136-0. IkeV 
j ion 

^0.0136-13.6keV 


emission in the bolometric luminosity LboI = lO^^ " erg s“' for 
NGC 3227. Their lower black hole mass results in a higher Ed¬ 
dington ratio of Tedd = 0.039. We find via linear interpolation, 
0.23 LboI ~ 10 "'^-® ergs^', which results in 
1x10''^ erg s“' using our measured X-ray data. 
Considering the uncertainties of this method, this value is con¬ 
sistent with our result from a more complete SED in a single 
source study. 


5.3.1. Constraints from the ionization parameter 

We assume that the illuminated cloud is in photoionization equi¬ 
librium with the radiation field. Using the definition of ^ (Eq. 1) 
and estimating the particle density from the radial extent, AR 
(< R), of the cloud, hh = Nh/AR, we find 


R < 




(3) 


This yields an upper limit for the distance of the cloud of 2.2 pc 
(2938 Id) when using the average measured parameters. 

5.3.2. Constraints from the photo-ionization equilibrium 


Following [Reynolds & Fabian| ( |1995] l, we now combine the in¬ 
formation we gain from Eq. 3 with requirements for the recom¬ 
bination time scale. As we assume that the absorber is in photo¬ 
ionization equilibrium with the radiation field, the recombination 
timescale of the plasma must be much smaller than the timescale 
of variations of the ionizing continuum. A direct test of this as¬ 
sumption would be measurements of a positive correlation be¬ 
tween log^ and Lion or Lx. Figure [T4| indeed shows a positive 
trend, but the correlation is weak. This can be explained with the 
assumption of a spherical cloud, which may be an oversimpli¬ 
fication. Any variations of the extend or number density of the 
cloud can distort a direct correlation. Also, the ionizing luminos¬ 
ity can not be measured directly, as the gap between the far-UV 
and soft X-rays is not observable. If the cloud reaches ionization 


equilibrium, the recombination time-scale is given by (Blustin 
|etal.|2005l[Krolik|19991 l 


(rieCki-ec) 


2 X 


(4) 


where ng = riellff cm ^ with «e being the electron particle den¬ 
sity, and where a^ec is the recombination rate coefficient, Z the 
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atomic number of the ion, and - T/IO^K ~ 1 , a represen¬ 
tative value corresponding to gas with the ionization parame¬ 
ter log^ > 1 ( Reynolds & Fabian||1995 Krolik & Kriss|[20OT] 


[Chakra vorty et al.||2009| l. We use Z = 9 to represent the likely 

dominant ions for WA 2 , Oviii and Neix (|Ka Iman & Bautista 


2001 1 . In their study of warm absorber properties, Chakravorty 


et al. 


(2009 1 find that for typical AGN SEDs the log^/T-logT 


stability curves are independent of the hydrogen number den¬ 
sity for 10^ < tin < 10'°cm“^. Assuming solar abundances, 
He = 1.2nH- Using these assumptions we find free ~ 370rig^ s, 
suggesting recombination timescales easily less than hours-days 
for most densities of relevance here. 

A further constraint on the density can be obtained from the 
observed variability of the absorber. The shortest timescales over 
which we can reliably measure changes in the illuminating flux 
and ionization in NGC 3227 is 7 days. RXTE light-curves indi¬ 
cate that NGC 3227’s X-ray flux typically varies by factors of 
a few tens of percent on timescales of 7 days and less ( jUttley 
|& McHardy||2005] l, which makes the assumption of photoion¬ 
ization balance a reasonable assumption. The weekly spaced 
Suzaku observations confiiTn this finding, while no variability 
is found within each of the observations except of Suz 1. The 
overall low variability is consistent with the assumption of photo 
ionization balance, such that free < Em - 7 d. This requirement 
can be translated to find a lower-limit on the hydrogen number 
density from Eq. 4, 


> 


2x 10^ 

1.2Z2fvar 


1.94 X 10^ cm^^ 


(5) 


which, in turn, sets an upper limit on the thickness of the absorb¬ 
ing layer. 


^R < —. 


( 6 ) 


EigurefTSjhelps to find an appropriate upper limit on the distance 
R of the absorber by considering the variation of log AR against 
log/? (see also jReynolds & Eabian|1995]l. We can exclude three 
regions in the log/?-log A/? space (Eig. HD light-red region). The 
first region is obtained by using Eq. 6 with the assumed average 
Ah of the absorbing cloud and the minimum hydrogen number 
density necessary to obtain ionization balance. In addition we 
can also exclude a region where the basic assumption A/? < /? is 
violated. Since 


^R = 


Ah/?"^ 


(7) 


and because of Ah = uh A/? we can also find an upper limit for R. 
Inserting the average Ah, ionization state, and luminosity mea¬ 
sured with Suzaku gives the dashed line in Eig. 15 An upper 


limit for R can be derived from the intersecting point of Eq. 6 
and Eq. 7. As a result, the blank region includes all values of AR 
and R that come into consideration based on our spectral anal¬ 
ysis, i.e., R < 10^*'^ cm = 0.65 pc = 770 Id (blue upper limit 
in Eig. Ml consistent with the dusty torus following jSuganuma] 
et al. ( 12006 1 , Ramos Almeida et al. ( 201 Ij ) and our geometrical 
considerations in Eq. 3. 


5.3.3. Constraints from a putative orbiting cloud 

In this section we additionally consider the Keplerian orbital 
motion of an obscuring cloud that is illuminated by the central 
source while passing the line of sight in >35 d, equal to the dura¬ 
tion of the campaign. We again use the average parameters Ah, 

Article number, page 12 of 



logR [cm] 

Fig. 15. Constraints on the radial distance R of the ionized absorbing 
material from the central engine visualized in the log R - log AR plane, 
where AR resembles the diameter of the intrinsic absorber. The regions 
filled with light-red can be excluded based on spectral results of Suzaku 
and the relations denoted along the border lines. The dashed line follows 
the relation A/? = AH/?^^/Lioii with the appropriate average parameters 
derived from spectral fits to all Suzaku observations. The intersecting 
points with line where AR < Nyilnn marks the upper limit on R. 


log^, and Lion. With these assumptions we can estimate the dis¬ 
tance of such a cloud. 

The first rough estimate is based on [Risaliti et al. ( 2007] 
|2009a|bj l and jLohfink et al.| ( |2012| l, who discuss a spherical cloud 
that passes the line of sight. It is able to fully cover the central 
X-ray source, if it is larger than the diameter of the X-ray emit¬ 
ting corona ~10rg ([Morgan et al.|2012|l to five times this value, 
which is arbitrarily chosen. This results in a lower limit for the 
distance of such a cloud of R > 0.1 pc, marked as green lower 
limit in Eig. [^ Due to the complex shape of the Ah profile, 
the underlying assumptions are probably too simple. The limit is 
also highly sensitive to the choice of the size of the X-ray emit¬ 
ting region. 

This very rough estimate of the distance can be signifi¬ 
cantly improved when adding the information on the average 
irradiating luminosity as well as column density and ioniza¬ 
tion state of the putative cloud. Writing the Kepler velocity as 
V - ^JG MbhIR - ARjAt - Ah/ (ny^Af) and expressing hh in 
terms of the ionization parameter gives ([Lamer et al.|2003[). 


/? >4x lO^^cmM^^ 


Tion ^Uays \ 

N22^ I 


2/5 


( 8 ) 


With M7 = Mbh/IO'Mq = 0.76, L 42 = L/10^'^ergs“‘ = 8.9, 
A 22 = An/lO^^cm^^ = 11.18, and At > 35 d we find R > 
0.017 pc = 20.7 Id. This lower limit is shown as additional green 


lower limit in Eig. 16 


Einally, we estimate the size AR of the assumed spherical 
cloud. Equating the Keplerian velocity with ARjAt = Nufna At, 
we find 


AfVGMBH' 


(9) 


Our distance estimates of 0.017-0.65 pc and the campaign- 
average value for Ah, then yield number densities of tin ~ 
2.7 X 10^-1.7 X 10^ cm~^, consistent with the lower-limit found 
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from the recombination time-scale analysis. This result trans¬ 
lates to a range in the size of the absorbers of AR ~ 6.6 x 10^^- 
4.1 X10*"* cm. Note that the column density profile has no defined 
ingress or egress. If we assume that the absorption by the cloud 
takes < 2 years, i.e., the time interval between the absorbed ob¬ 
servations and the 2006 unabsorbed XMM-Newton observation, 
the density would be about one order of magnitude less and its 
size accordingly larger. 

5.4. A dust-free cloud located In the BLR? 

The range of 0.017-0.65 pc found in the previous section means 
the possible location for the cloud spans radial distances both 
inside and outside the dust sublimation zone. Here, we define 
the dust sublimation zone to range from 0.4Rd-f^d, with the the 
dust sublimation radius. 


/ Lboi Td ® 2.5 X 10**-1.5 X 10"* cm is consistent with the typical dens: 

Rd - 0.4pc ■ I j ^45 I ( ~ P** >10^ cm~^ expected for BLR clouds (Baldwin et al.||2003 


due to a X-ray absorbing cloud that fully covers the X-ray corona 
but not the entire UV-emitting disk. 

A second possible reason for the low reddening is that the 
cloud does indeed contain dust but at an extremely high gas-to- 
dust ratio. With the results from above, we find a gas-to-dust 
ratio of AhMv ~ 2.5 x 10^^ cm~^ mag^*, which is about a factor 
100 higher than the assumed Galactic gas-to-dust ratio. 

Finally, the low level of variability that we infer for the red¬ 
dening supports the notion that the reddening matter has an ori¬ 
gin that is distinct from the much more strongly variable X-ray 
column. It is consistent with the distant, dusty, “lukewarm” ab¬ 
sorber identihed by [Kraemer et al.| ( |2000) l and [Crenshaw et al.| 
pOOlj l as well as recent Swift observations. They therefore con¬ 
clude that the X-ray absorbing gas is likely dust-free, which is 
supported by Sect. 5.3.3| We find that the absorber can be lo¬ 
cated well below the dust sublimation radius, and thus plau¬ 
sibly be a BLR-cloud. The hydrogen number density range of 
2.5 X 10^-1.5 X 10® cm“^ is consistent with the typical density of 


Os- 


143 c 1 and 


with an assumed bolometric luminosity of 1.3 x lO"'" erg s 
a dust temperature of Td = 1500 K (|Barvainis||1987| [Nenkova 
jet al.|2008b| l. _ _ 

The result is consistent with [Blustin et al.j P005| l, who de¬ 
termine the minimum and maximum distance of the absorbing 
gas solely based on geometrical considerations for a sample of 
23 Seyfert galaxies of intermediate classification. [Blustin et al.j 
find that most low-velocity absorbers are consistent with the in¬ 
ner edge of the torus. Similarly, [Risaliti et al. ( 2002[ ) explain the 
column density variations of a large fraction of their Seyfert 2 
sample sources with clouds at the inner edge of the dust subli¬ 
mation zone. 

We can test the dust content of the cloud by using the redden¬ 


ing derived in Sect. 4.1 If the Galactic dust-to-gas ratio is appli¬ 


cable, the X-ray absorbing gas columns predict that we should 
see roughly Ay ~ 20-60 mag of optical extinction. During the 
strong 2008 absorption event, however, we only measure a red¬ 
dening of Ay ~ 0.45 mag ~ 1.2 x 10^' cm“^. 

There are several potential reasons for the lack of a strong 
dust component in the variable X-ray column density. First we 
consider a scenario where a cloud inside the line of sight to the 
X-ray source may indeed contain dust, but does not cover the line 
of sight to the UV continuum source. The diameter of a spherical 
cloud on a Keplerian orbit was determined in Sect. |5.3.3[ to be 
~6.6x lO'M.l X 10'"* cm. We now estimate the diameter of the 
UV emitting part of the accretion disk for comparison. The radial 
temperature profile of a standard thin disk is 


T(R) ~ 2 X 10^ K 


M 


lO^M© 


\MEdd/ U0*"*cm/ 


with the derived Eddington ratio Tedd = i'boi/^'Edd = MjM^di ~ 
0.004, the black hole mass and the radius of the UV emitting re¬ 
gion. To solve the equation for R, we estimate the temperature 
at the outer UV emitting disk using the Wien displacement law 
and the longest wavelength Rwi/f/UVOT filter. We find a radius 
of 1.1 X 10*^ cm, i.e., larger than the estimated diameter range 
of a spherical, homogeneous X-ray absorbing cloud. In addi¬ 
tion, McHardy et al.[ ( 2014| l, studying the interband time lags in 
NGC 5548, inferred that UV-emitting regions can extend slightly 
further than expected from standard cr-disk theory, at least in that 
object. Consequently, based on our results that assume a simple 
spherical cloud, we cannot rule out that the weak reddening is 


[terbrock|1989| l. 

5.5. The morphology and spatial distribution of a putative 
cloudy absorber 

We now investigate the possibility for the detected absorber to 
be part of an overall clumpy medium. 


According to the clumpy torus model of Nenkova et al. 
P008a|b[ ), the average number clouds along a line of sight with 
inclination i with respect to the system symmetry axis, Nc, is 
given by 


Nc = No exp 


90-i 


■\2 


( 12 ) 


where No is the number of clouds along a ray in the equato¬ 
rial plane between the Rd and YRd. The radial distribution of 
the clouds follows a power law This model has observa¬ 
tional support from extensive time-resolved X-ray spectroscopy 
(Markowitz et al. 2014|l an d Bayesian model hts (Asensio Ramos 
& Ramos Almeida 2()09|l to IR SEDs ( [Ramos Almeida et al. 
201 l[|Alonso-Herrero et al.|201 1| |. 


Elitzur ( 2001) and Gaskell et al.[P008l l claim that clouds can 
exist both in the BLR below, and in the dusty torus beyond Rd, all 
as part of a common structure. We therefore assume for simplic¬ 
ity that the power-law index q is the same for both the BLR and 
the torus clouds, which yields a successively growing number of 
clouds towards the center. We assume the following mode values 
for the clumpy torus parameters from the posterior probability 
distributions of |Alonso-Herrero et al.[ ( |T011| l: Nq = 15, T = 15, 
q - O.l, (T - 44° and i - 30°. They use IR photometry and ad¬ 
ditionally spectroscopic data around 10//m. These values result 
in ~16 clouds in the equatorial plane between 0.4 and 15Rd. The 
lower limit equals the derived minimal distance of the absorber 
in NGC 3227. We call this number V^^***, as it includes both 
dusty clouds detected in the IR and also dust-free clouds that ad¬ 
ditionally absorb the X-rays. The number of clouds at the given 
inclination angle and width of the Gaussian cloud distribution 
then equals = 2.5 (Eq. 12). 

Eor comparison one can estimate the average number of 
clouds on the line of sight based on 6.9 years of sustained 
RXTE monitoring, the 2006 XMM-Newton observation, and the 
35 days of Suzaku and Swift monitoring in 2008. During that 
time NGC 3227 was in an obscured state for a total of at least 
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Fig. 16. Radial distribution of clumps that can potentially exist below and beyond the dust sublimation zone, i.e., (dashed lines). The 

distance limits from Sect. |5.3| are marked with colored arrows. For all estimates we use the average measured parameters of WA 2 plus the average 
ionizing luminosity. The red upper limit from Sect. |5.3.T] uses the definition of the ionization parameter, while the blue upper limit from Sect. 5.3.2| 
is determined from geometrical constraints on a spherical homogeneous cloud in photoionization balance. The Keplerian orbital motion of such 
a hypothetical cloud passing the line of sight is additionally included to form the lower-limit in green (Sect. |5.3.3^ while the second, larger 
lower-limit in green is determined for an orbiting cloud neglecting the information on its ionization or column density. 


~114 days. Based on the given data, NGC 3227 spends 4.4% of 
the observed time in eclipse. 

Assuming Poisson statistics, the time averaged escape prob¬ 
ability for photons without undergoing strong absorption by an 
average of clouds in the line of sight is given by 


P ~ exp(-A^+'’^T.r) 


(13) 


where the monochromatic optical depth of one cloud, > I 
in order to be able to obscure the line of sight as observed dur¬ 
ing eclipses. With a 4.4% probability of obscuration, the inte¬ 
grated escape probability equals ~95.5%, i.e., ~ 0.045. 

This value is clearly lower than the average of 2.5 clouds esti¬ 
mated above. 


One can still find consistency when considering severe de¬ 


generacies between different attempts of clumpy tori fits. Ramos 
Almeida et al. (201 l|l for example only use IR photometry be 


tween 1-40/tm. They find much larger uncertainties and differ¬ 
ent mode values of Aq = 2, T = 19, q - 0.6, a - 33° and 
i = 66°. Based on these values we find = 2.3 between 0.4 
to 30Rd; is consequently equal to 1.4. When inserting not 

the mode values but other values of high probability within the 
posterior probability distributions of the parameters, we can find 
consistency with the observed number of 0.045 clouds along the 
line of sight. If the inclination angle is even more face-on, the 
number of clou ds A^^ ^^ on the line of sight is also significantly 
reduced.[Marin|(|^14[) provides an overview over different meth 


ods and inclination angles measured so far. Hicks & MaUcan 
( |2008| l and [Fischer et al.|P013| l indeed find a value of i ~ 16° 
using near-IR spectroscopy consistent with [Fischer et al.|(|2013| l 
although the results from the optical polarimetry (Sect. [4.2| l are 
more consistent with the intermediate inclinations derived from 
IR SED fitting. 


We note that the fact that we obtain a larger number of clouds 
in the line of sight from extending the clumpy torus model than 
required from historical eclipse events is not inconsistent with 
the spectral results. [Markowitz et al.[ ( [2009| l show that even in 
a relatively unabsorbed state, ionized absorbers are present in 
the line of sight. Those could explain the determined excess of 
clouds. 
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5.6. The origin of the absorber 

Here, we discuss the implications of the measured variable den¬ 
sity profile of the absorber. Besides the strongly variable column 
density over time, the covering fraction remains roughly constant 
for the latter five out of six Suzaku observations. This argues 
against two distinct clouds as explanation for the two overden¬ 
sities but for a potential filament with a variable internal density 
structure that enters the line of sight. As a test, we can estimate 
the limiting size of one cloud, which is confined by its own mag¬ 
netic field or the external pressure of the intercloud medium, not 
to get tidally sheared by the central force of gravity, to 


rci < 10'® 


M^ 


cm 


(14) 


(Elitzur & Shlosman||2006 1 where A 23 = AH/10^^cm“^, Rp^ - 
jf/lpc IS the distance of the cloud to the central engine, and 
where Mj = Mbh/IO’Mq. Inserting the average column den¬ 
sity Ah from the 2008 event and Rp^. ~ 0.017-0.65, we find 
rd < 3.3 X lO*® - 1.8 X 10'®cm. Eor comparison, the esti¬ 
mated diameter of a spherical cloud on a Keplerian orbit was 
6.6 X lO'^cm^.l X 10 '^ cm. Both ranges are consistent. Note, 
however, that the limiting size of the absorber, r^, is strongly 
dependent on its distance. We can argue for a distance Rp^ of 
the absorbing cloud towards the lower limit of allowed ranges 
when considering the weak reddened, implying a location inside 
the dust sublimation radius. In this case we would have to as¬ 
sume a rather extended and filamentary medium in contrast to 
the simplified picture of spherically symmetric clouds bound by 
self-gravity. 

Mechanisms that are thought to be able to lead to such struc¬ 
tures include MHD-driven winds that are launched from the ac¬ 
cretion disk, for example, via magneto-centrifugally accelera¬ 


tion ( Blandford & Payne|1982[ Contopoulos & Lovelace|1994 

[Konigl & Kartje[[1994| l. We find consistency between the mea¬ 

sured ionization structure and column density of the absorber 
and the range of values simulated for MHD winds by[Fukumura 
leTaLldT^ . 


The absorber is also qualitatively consistent with the turbu¬ 
lent, dusty disk wind proposed by [Czerny & Hryniewicz[ ( |2011[ l. 
In this model, dusty clouds from the low-ionization part of the 
BLR can rise from the disk, e.g., driven by MHD-forces, and 
get irradiated and heated sufficiently such that dust sublimates. 
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Eventually clouds may fall back towards the disk when the grav¬ 
itational force prevails over the driving forces. In this context the 
absorbing cloud from 2008 may be just recently sublimated and 
temporarily situated at a height intercepting our line of sight. 

The question remains how such driven clouds can get 
stretched out to filamentary structures of internal overdensities 
that seem to best explain the observations. A potential mecha¬ 
nism is described by [Emmering et al.| ( |1992| l in the context of 
MHD winds. As a wind or cloud forms, the material initially 
contains a high amount of dust, which is a good coolant via IR 
thermal emission. Once the clouds are driven away from the ac¬ 
cretion disk they are heating up to ~10^ K as the dust sublimes. 
After a cloud is photo-ionized, its internal pressure increases 
and expands at roughly its internal sound speed until its pres¬ 
sure decreases to roughly that exerted by the external magnetic 
field. Rees ( |1987| l show that the polodial magnetic field of the 
disk together with the internal field of magnetized BLR clouds 
can cause the expanded clouds to stretch out along the magnetic 
field, which is frozen in the disk. 

Additionally we mention radiative driving of accretion disk 
winds. Dorodnitsyn & Kallman ( 2012| l provide detailed simu¬ 
lations on how the X-ray radiation from the hot corona can be 
transferred to the IR by exerting pressure on dust and allowing 
clumpy or filamentary structures to arise. Again, the cloud or fil¬ 
ament can be ionized by the nuclear radiation when reaching a 
certain height matching the observables of the absorber. 

The 2008 absorption event might not be unique. Similar to 
the long-term obscurations of ~80d and > 35d of NGC 3227 
in 2000/2001 and 2008, respectively, NGC 5548 is continuously 
obscured for a few years ( [Kaastra et al.|2014) l. Such phenomena 
may thus be typical to Seyfert galaxies although the probability 
to observe them is low. 


6. Summary and Conclusions 

We have studied data from a 5 week long Suzaku and Swift X- 
ray and UV monitoring campaign on the Sy 1.5 NGC 3227 in 
late 2008, which caught the source in an absorbed state. We per¬ 
formed a time-resolved X-ray spectroscopy to untangle the var¬ 
ious emission components, i.e., the coronal power law and ion¬ 
ized reflection, as well as multiple ionized absorption compo¬ 
nents, and isolate the time-dependent behavior of the dominant 
absorber in 2008. Past RXTE monitoring reveals NGC 3227 to 
be usually unabsorbed, except for eclipses in 2000/2001 (Lamer 
|et al.|20()^ and 2002 ( [Markowitz et al.|2014| l. The ultimate aim 
of this study is to understand the properties of the X-ray ab¬ 
sorber. 

We used two archival and previously published XMM- 
Newton spectra to derive a baseline model that we fit simultane¬ 
ously to six Suzaku and twelve Swift observations from 2008. An 
extensive exploration of the parameter space reveals only minor 
degeneracies. The two XMM-Newton spectra allowed us to iden¬ 
tify three differently ionized absorbers, called WAj (lowly ion¬ 
ized), WA 2 (intermediately ionized) and WA 3 (highly ionized). 
We conclude that the detected absorption variability in 2008 is 
solely due to WA 2 , while WAl and WA 3 remain constant. 

We resolve the column density and covering fraction of WA 2 
to be variable with time. The covering fraction rises from ~70% 
for the first Suzaku observation up to ~90% for the remain¬ 
ing five observations. The column density varies from ~5 to 


-18 X 10^^ cm“^ in a doubly-peaked density profile. The ioniz¬ 
ing luminosity Lion and the ionization state ^ of WA 2 are seen to 
vary by factors of 1.3 and 104, respectively. A correlation anal¬ 


ysis shows a tentative dependence between f and Lion as a proxy 
for Lion. 

We use the average parameters of WA 2 ((Nh) ~ 11.2 x 
lO^^cm^^, (log^) ~ 1.1, (Lion) ~ 8.9 x 10"^^ergs“') during the 
>35 d long absorption event to estimate that the distance of the 
absorbing medium is between 0.017 and 0.65 pc from the central 
SMBH. Here we made use of all available information, which 
comprise the ionization state, column density and minimal du¬ 
ration of the event as well as the incident irradiation. The un¬ 
derlying, simplified assumption is that of a spherical cloud on a 
Keplerian orbit around the illuminating source. 

Eor the derived distance range we infer the hydrogen number 
density of the spherical cloud to be ~ 2.7 x 10^ - 1.7 x 10® cm“^. 
The cloud therefore has a diameter of AR ~ 6.6 x 10*^ - 4.1 x 
10 '"^ cm. We conclude that the absorber may be located in the 
outermost dust free BLR or the inner dusty torus. 

We And only moderate reddening in the Swift UVOT data, 
which is due to a column density of Ah ~ 1.2x10^^ cm“^. In con¬ 
trast, we measure X-ray-absorbing column densities being about 
a factor of 100 larger. We conclude based on Sect.j^that the X- 
ray absorber responsible for the variability in 2008 and the ab¬ 
sorber responsible for extinction are distinct. The first seems to 
comprise a high gas-to-dust ratio due to its vicinity to the SMBH. 
A location in the BLR would be consistent with our results. The 
latter is consistent with the distant (100 pc) lukewarm absorber 
proposed by Kraemer et al. ( 2000| ) and Crenshaw et al. ( 2001| l, 
which still seems to be present in recent Swift observations. 

We also investigated the measured absorber in the context 
of the clumpy absorber model of Nenkova et al. ( 2008a|b l. The 
inferred distribution of clouds is consistent with the eclipses ob¬ 
served for NGC 3227. The result that the absorber is likely lo¬ 
cated in the BLR, leads us to extend the distribution of clouds 
also down to 0.4Rd- The consistency of the model with past ob¬ 
servations and the 2008 event is in particular reached for small 
inclination angles (< 20-30°) of the putative clumpy torus that 
are measured by independent methods and help to explain the 
predominant lack of absorption in NGC 3227 during past obser¬ 
vations. 

A spherical cloud situated below the dust sublimation zone 
can most likely not withstand tidal shearing by the central source 
of gravity. Several mechanisms including MHD-driving are able 
to lift an initially dusty cloud from the disk where it is getting 
exposed to nuclear radiation and therefore ionized. Pressure gra¬ 
dients potentially stretch the cloud preferably along the magnetic 
held lines frozen to the disk. 

In that context, the observing campaign represents a series of 
snapshots of parts of a cloud or filament whose properties (size, 
density distribution) may be evolving over time. The density 
profile of this event contrasts with that of the 2000/2001 event, 
which showed a centrally symmetric profile and clear ingress 
and egress. We therefore conclude that NGC 3227 is a rare lab¬ 
oratory to study the range of physical processes that form and 
shape clouds or filaments. 

Acknowledgements. This research was partially funded by a grant from the Bun- 
desministerium fur Wirtschaft und Technologie under Deutsches Zentrum fur 
Luft- und Raumfahrt grant 50 OR 1311. We thank Mirko Krumpe and Robert 
Nikutta for reading the manuscript and giving very helpful comments. We made 
use of ISIS functions provided by ECAP/'Remeis observatory and MIT i http: 
//www.sternwarte.uni-erlangen.de/isis/i as well as the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet Propulsion Labora- 
toiy, California Institute of Technology, under contract with the National Aero¬ 
nautics and Space Administration. This work used data obtained with the Suzaku 
satellite, a collaborative mission between the space agencies of Japan (JAXA) 
and the USA (NASA), as well as data from the Swift satellite, a NASA mission 
managed by the Goddard Space Flight Center, XMM-Newton, an ESA science 
mission with instruments and contributions directly funded by ESA Member 


Article number, page 15 of 19 












































A&A proofs: manuscript no. ngc3227 


States and NASA, the Wide-held Infrared Survey Explorer, which is a joint 
project of the University of California, Los Angeles, and the Jet Propulsion 
Laboratory/California Institute of Technology, funded by NASA, and the Two 
Micron All Sky Survey, a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/Califomia Institute of Technology, 
funded by NASA and the National Science Foundation. We also used observa¬ 
tions made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Science Institute. STScI is operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. The research is also based upon observations acquired at the Ob- 
servatorio Astronomico Nacional in the Sierra San Pedro Martir (OAN-SPM), 
Baja California, Mexico. TM acknowledges support by UNAM-DGAPA Grant 
PAPIIT IN104113 and CONACyT Grant Cientfhca Basica 179662. We thank J. 
E. Davis for the development of the slxfig module that has been used to prepare 
the figures in this work. 


References 

Akylas A., Georgantopoulos I., Griffiths R.G., et al., 2002, MNRAS 332, L23 

Alonso-Herrero A., Ramos Almeida C., Mason R., et al., 2011, ApJ 736, 82 

Antonucci R., 1993, ARA&A 31, 473 

Arav N., Barlow T.A., Laor A., et al., 1998, MNRAS 297, 990 

Asensio Ramos A., Ramos Almeida C., 2009, ApJ 696, 2075 

Baldwin J.A., Ferland G.J., Korista K.T., et al., 2003, ApJ 582, 590 

Barvainis R., 1987, ApJ 320, 537 

Bentz M.C., Peterson B.M., Netzer H., et al., 2009, ApJ 697, 160 

Bessell M.S., 2005, ARA&A 43, 293 

Bessell M.S., Castelli R, Plez B., 1998, A&A 333, 231 

Beuchert T., Wilms J., Kadler M., et al., 2013, ArXiv e-prints 

Blandford R.D., Payne D.G., 1982, MNRAS 199, 883 

Blustin A.J., Page M.J., Fuerst S.V., et al., 2005, A&A 431, 111 

Bun'ows D.N., Flill J.E., Nousek J.A., et al., 2005, Space Sci. Rev. 120, 165 

Calderone G., Sbarrato T., Ghisellini G., 2012, MNRAS 425, L41 

Chakravorty S., Kembhavi A.K., Elvis M., et al., 2009, MNRAS 393, 83 

Contopoulos J., Lovelace R.V.E., 1994, ApJ 429, 139 

Crenshaw D.M., Kraemer S.B., Bruhweiler F.C., Ruiz J.R., 2001, ApJ 555, 633 
Czerny B., Hryniewicz K., 2011, A&A 525, L8 
Davies R.I., Thomas J., Genzel R., et al., 2006, ApJ 646, 754 
de Vaucouleurs G., de Vaucouleurs A., Corwin, Jr. H.G., et al., 1991, Sky Telesc.. 
82, 621 

Dorodnitsyn A., Kallman T., 2012, ApJ 761, 70 

Elitzur M., 2007, In: Ho L.C., Wang J.M. (eds.) The Central Engine of Active 
Galactic Nuclei. Astron. Soc. Pacihc Conf. Ser. 373, Astron. Soc. Pacihc, San 
Francisco, p. 415 

Elitzur M., Shlosman I., 2006, ApJ 648, LlOl 

Elvis M., Risaliti G., Nicastro R, et al., 2004, ApJ 615, L25 

Emmering R.T., Blandford R.D., Shlosman I., 1992, ApJ 385, 460 

Eischer T.C., Crenshaw D.M., Kraemer S.B., Schmitt H.R., 2013, ApJS 209, 1 

FitzpatrickE.L., 1999, PASP 111, 63 

Fukumura K., Kazanas D., Contopoulos I., Behar E., 2010, ApJ 715, 636 
Garcia J., Dauser T., Lohhnk A., et al., 2014, ApJ 782, 76 
Garcia J., Dauser T., Reynolds C.S., et al., 2013, ApJ 768, 146 
Gaskell C.M., Goosmann R.W., Klimek E.S., 2008, Mem. Soc. Astron. Italiana 
79, 1090 

George I.M., Turner T.J., Netzer H., et al., 1998, ApJS 114, 73 
Ghisellini G., Maraschi L., Tavecchio R, 2009, MNRAS 396, L105 
Goiford J., Reeves J.N., Tombesi R, et al., 2013, MNRAS 430, 60 
Guainazzi M., Nicastro R, Fiore R, et al., 1998, MNRAS 301, LI 
Heman-Caballero A., Alonso-Herrero A., Hatziminaoglou E., et al., 2015, ArXiv 
e-prints 

Hicks E.K.S., Malkan M.A., 2008, ApJS 174, 31 

Houck J.C., Denicola L.A., 2000, In: Manset N., Veillet C., Crabtree D. (eds.) 
Astronomical Data Analysis Software and Systems IX. Astron. Soc. Pacihc 
Conf. Ser. 216, Astron. Soc. Pacihc, San Erancisco, p. 591 
Kaastra J.S., Kriss G.A., Cappi M., et al., 2014, ArXiv e-prints 
Kalberla P.M.W., Burton W.B., Hartmann D., et al., 2005, A&A 440, 775 
Kallman T., Bautista M., 2001, ApJS 133, 221 
Kinney A.L., Calzetti D., Bohlin R.C., et al., 1996, ApJ 467, 38 
Komossa S., Pink H., 1997, A&A 327, 483 
Konigl A., Kartje J.R, 1994, ApJ 434, 446 
Koratkar A., Blaes 0., 1999, PASP 111,1 
Koyama K., Tsunemi H., Dotani T, et al., 2007, PASJ 59, 23 
Kraemer S.B., George I.M., Turner T.J., Crenshaw D.M., 2000, ApJ 535, 53 
Krolik J.H., 1999, Active Galactic Nuclei, Princeton Univ. Press, Princeton, NJ 
Krolik J.H., Kriss G.A., 2001, ApJ 561, 684 
Ki'ongold Y., Nicastro P., Elvis M., et al., 2007, ApJ 659, 1022 
Kuhnel M., Muller S., Kreykenbohm L, et al., 2015, Acta Pol. 55, 1 
Lamer G., Uttley P, McHardy I.M., 2003, MNRAS 342, L41 


Article number, page 16 of 


Landt H., Bentz M.C., Ward M.J., et al., 2008, ApJS 174, 282 
Lohhnk A.M., Reynolds C.S., Mushotzky R.R, Wilms J., 2012, ApJ 749, L31 
Maiolino R., Risaliti G., Salvati M., et al., 2010, A&A 517, A47 
Marin R, 2014, ArXiv e-prints 

Markowitz A., Reeves J.N., George I.M., et al., 2009, ApJ 691, 922 
Markowitz A.G., Ki-umpe M., Nikutta R., 2014, MNRAS 439, 1403 
McHardy I.M., Cameron D.T., Dwelly T., et al., 2014, MNRAS 444, 1469 
McKeman B., Yaqoob T., 1998, ApJ 501, L29 
Morgan C.W., Hainline L.J., Chen B., et al., 2012, ApJ 756, 52 
Mould J.R., Huchra J.P., Freedman W.L., et al., 2000, ApJ 529, 786 
Munoz Marin V.M., Gonzalez Delgado R.M., Schmitt H.R., et al., 2007, AJ 134, 
648 

Mundell C.G., Holloway A.J., Pedlar- A., et al., 1995a, MNRAS 275, 67 
Mundell C.G., Pedlar A., Axon D.J., et al., 1995b, MNRAS 277, 641 
Nenkova M., Ivezic Z., Elitzur M., 2002, ApJ 570, L9 
Nenkova M., Sirocky M.M., Ivezic Z., Elitzur M., 2008a, ApJ 685, 147 
Nenkova M., Sirocky M.M., Nikutta R., et al., 2008b, ApJ 685, 160 
Nowak M.A., Neilsen J., Markoff S.B., et al., 2012, ApJ 759, 95 
Osterbrock D.E., 1989, Astrophysics of gaseous nebulae and active galactic 
nuclei. University Science Books, Mill Valley, CA 
Peterson B.M., Ferrarese L., Gilbert K.M., et al., 2004, ApJ 613, 682 
Quillen A.C., Alonso-Hen'ero A., Rieke M.J., et al., 1999, ApJ 527, 696 
Ramos Almeida C., Levenson N.A., Alonso-Herrero A., et al., 2011, ApJ 731, 
92 

Rees M.J., 1987, MNRAS 228, 47P 
Reynolds C.S., Fabian A.C., 1995, MNRAS 273, 1176 
Risaliti G., Elvis M., Fabbiano G., et al., 2007, ApJ 659, LI 11 
Risaliti G., Elvis M., Nicastro R, 2002, ApJ 571, 234 
Risaliti G., Miniutti G., Elvis M., et al., 2009a, ApJ 696, 160 
Risaliti G., Nardini E., Salvati M., et al., 2011, MNRAS 410, 1027 
Risaliti G., Salvati M., Elvis M., et al., 2009b, MNRAS 393, LI 
Rivers E., Markowitz A., Rothschild R., 201 la, ApJ 742, L29 
Rivers E., Markowitz A., Rothschild R., 2011b, ApJS 193, 3 
Roming P.W.A., Kennedy T.E., Mason K.O., et al., 2005, Space Sci. Rev. 120, 
95 

Skrutskie M.R, Cutri R.M., Stiening R., et al., 2006, AJ 131, 1163 

Smith J.E., Robinson A., Alexander D.M., et al., 2008, MNRAS 350, 140 

Stalevski M., Fritz J., Baes M., et al., 2012, MNRAS 420, 2756 

Strirder L., Briel U., Dennerl K., et al., 2001, A&A 365, L18 

Suganuma M., Yoshii Y, Kobayashi Y, et al., 2006, ApJ 639, 46 

Takahashi T, Abe K., Endo M., et al., 2007, PASJ 59, 35 

Tarter C.B., Tucker W.H., Salpeter E.E., 1969, ApJ 156, 943 

Tombesi R, Cappi M., Reeves J.N., et al., 2013, MNRAS 430, 1102 

Turner M.J.L., Abbey A., Arnaud M., et al., 2001, A&A 365, L27 

Turner T.J., Reeves J.N., Kraemer S.B., Miller L., 2008, A&A 483, 161 

Uttley R, McHardy I.M., 2005, MNRAS 363, 586 

Vasudevan R.V., Fabian A.C., 2009, MNRAS 392, 1124 

Vasudevan R.V., Fabian A.C., Gandhi P, et al., 2010, MNRAS 402, 1081 

Vemer D.A., Ferland G.J., Korista K.T., Yakovlev D.G., 1996, ApJ 465, 487 

Wilms J., Allen A., McCray R., 2000, ApJ 542, 914 

Winge C., Peterson B.M., Horne K., et al., 1995, ApJ 445, 680 

Wright E.L., Eisenhardt P.R.M., Mainzer A.K., et al., 2010, AJ 140, 1868 


19 



T. Beuchert et al.: A Variable-Density Absorption Event in NGC 3227 mapped with Suzaku and Swift 


» -a 

•“ I 

t/3 C 

8 I 

« s 

= I 

O ^ 
c 

C3 


s: 


cS o 


a '■ 

3 

cr 


ra = 
jq £ 


> 

tu 

3 
O 
o - 
3 ■ 

(D 

3 

o^ 7 

O ^ 

■ Oh 


C 

3 

.£ 

c 

tu 

> 

‘Sb 


33 D 
3 £ 

I ^ 

. Cl 
C 

N <D 

3 TD 
C/3 C 
— <u 

.2 <1^ 
■5 ’CJ 
> .£ 

^ OJ 

t/3 C3 

.3 
O C3 

«D TD 
3 

3 


>» cfl 
iZi r- 

^ O 

^ S 

O N 

c S 

s S 

o o 
* Z 


o 

o 

U4J 

73 


o 

cS 

73 

3 

72. 

'> 

'-B 


«D 


c 

3 


3 


3 

1°’ 
^ ffi ' 

° •§ - 
^ g 

3 M 
t/5 G . 

cfl 

C £Ji 

Q c 


3 . 

e 

«L) 

Cfl 

X) 

o 

O 

CX 

c/2 

3 

O 


£ X 

s s; 

P 


o 

S-H 

bfl , 
3 


3 3 
G cr 
.£ 


’B 2 ’ 

tH G ’ 

rP B 

U-i ■*-’ • 

C/2 

S3 ^ 
a3 

^ § 

3 3 
Oh CIh 

^ tS ' 
O C 

^ u 

c/2 

' 3 C 
^ D 
. Oh 
Tt <U 
0) ’? 

r3 £ 

H -x t 


/-N 


T—( 


'w/ 


0) 


1—1 


o. 


B 


■H 


{/) 

1 


1 


0) 


d 




+J 


* 


/-N 


(Nl 


'w/ 


pq 




u 


o. 


■H 


X 


N 


* 



NO 

»—( 

N 


3 

pq 

00 


u 


o. 


■H 


X 


N 


* 





in 


N 

+ 

3 

(Nl 

CO 

Vw/ 

pq 




u 


o. 


•H 


X 


N 

N 

* 

3 


CO 

Vw/ 


+J 


fC 

+J 


{/) 


d 


O 


U 

CO 

+ 

N 


3 


CO 

Vw/ 


+J 




fC 


+J 


(fi 


d 


O 


u 

CN 

1 

N 

1—( 

3 

Vw/ 

CO 

Vw/ 

* 








Vw/ 


U 


QJ 


> 


1—1 

r-H 

1—1 

N 

■H 

3 

X 

CO 

+ 





Vw/ 


1—1 


O. 




i+H 


O 


+J 




u 


Vw/ 


* 


/-N 


1—( 


Vw/ 


+J 


(fi 


d 


O 


U 


+J 


oi 


73 



time dependent 


— 00*0 mm o 

— OO O ■ ■ 


oodo—d'^ 3 -—'ood 
+ 1+ 1+ io-:+ 1+ I 
O CN O + I ^ ^ 

r-; r-; ^ ^ ^ 


ON 

o 

o 

,=== 

+1 

o 

m 

ON 


osmH~mfnr^r~ 
00003-——m 

dd+idddd 4-| 
+ I tJ- +1+1 
CM ON O 0^ 00 

00 NO cn NO r+ 


r- 
o 
o 
d 
+1 
ON 
<N 
(N ON 


S 

d 

+1 


■ cs 


*0 m 

WHm-^oom-^m asr^PP 
oood—odiN^m——OO 
+ idddddd—dd+ I 
IT) + i+ i+ i+ i+ [m 
^ <N ON ^ ON 00 m 
CO NO <N O O 00 


O (N 


o o 


NO 

o 

o 

d 


— mo 3" 
m cN f'' cs 

. .PP"rP3-oiO om 

pppppppcN^ +1 

d d 
+ i in 


o d + I + I d d 
+_ I CO nj +1 _|_| 


NO in ON CO 

r- NO ^ o <N 

d ^ d (N NO 


O \D 

CS bO O 


jodd^ 


C' CNl 

o — 
«o o 


+ I (N 
(N ON 

NO in (N 
d 


_OO O + I 

+ 1 + I + I CO 

3j- in o r- 


+1 O 0 - 

' I 4- I - 


O 

(N 


(N 00 

^ d 


O ?i 


— — O ^ 04 O 
_j_l I/-) — Tt rn 04 ^ 

d d d d d d +1 

00 + I + I + I 

in CO in ^ 

NO CO r- r- o 

—^ d ^ (N 


c/2 01 


3 " m 

04 04 

o o 
d d 


1 

o 

1 

o 

d 

o 

§ ^ 

X 

X 

£? 

M £0 

04 



o 

o sH- 

a 

a 


W ^ 

lu 

Li 

o 

c 

bJ) 

X bi) 

S: o 


d 




o 

d 


04 

PL 

d 

o 


S 




time-independent 


r- 

o 

o 

d 

+1 


NO 

o 

o 

d 

+1 


NO NO -f- 
ON CO NO 

00 ^ 


O 00 O — 
— OO — 

0000 
d dd d 


^ 

pp d 

o o 

+1 

SO ^ 

P d 
o I 


— —o o 
Tt fvj + I H- 

+ I r- o 

r-H p 

3t ^ r-H 


in 


00 



* 

r-> 

ON 

-t- 

OO m o 

O X + 1 

CO 


in 

in 

NO 

.—H 

CO 

CS 

o 

+ 1 ^ 

r- 

-1- 

q 

q 


00 

r-H 

d 

q 

P ^ 

d 

1-H 

1-H 



(N 

d 

1 


CO CO 


XXX 


Q 

X 

X 


N 


s 

o 




o 

£P 

<D 

bi) 

O 


o 

d 

o 




g + 




s 

O 

£? 

o 

•NNd 

Oi) 

O 




^ o 

p d 

NO 00 


o 

+1 

in 


o 

+1 
r- 

NO (N 


in 

o 

d 

+1 


in 
o 
d 

+1 CJ 

NO 00 


o 

d 

+1 


o 
d 

+1 

CO 
I> (N 


On 

O 


in 

o 

d 

+1 

o 

00 


CO 
CO O 

p d 
° +' 
+1 o 

NO NO 


!-h I 

+ ^ 
3 bD 
O ^ 

04 

>0 

I 

\o o 

2 7 

CO 'O 

2dXd 


X 

'■5 

c 

0) 

a 

a 


< 


bL 


Article number, page 17 of 


19 






















Table 5. List of parameters for the simultaneous fit of all Swift observations. The model equals the one fitted to the Suzaku data. The simultaneous fit consists of time-dependent parameters 
individually fitted to each observation and time-independent parameters that are equal for all observations, denoted with quotation marks. The symbol * marks parameters adopted from the 
simultaneous fit of the Suzaku observati ons. The redshift is fr ozen to the systemic value of z = 0.003859 (de Vaucouleurs et al.|1991). All spectra are additionally absorbed by a cold column of 
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T. Beuchert et al.: A Variable-Density Absorption Event in NGC 3227 mapped with Suzaku and Swift 
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Fig. 17. Spectra and model components of all Suzaku observations. We show the complete model (black solid line) with residuals, the hard X-ray 
power law (dark green dashed line), the ionized reflection component (red dotted line), the emission absorbed by the variable partially covering 
column (blue dotted-dashed line) and the leaked emission (cyan doubly-dotted-dashed line). 
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